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ON THE COMPLEX ROOTS OF 
ALGEBRAIC EQUATIONS* 


BY A. J. KEMPNER 


1. Introduction. The topic I have chosen is so vast that a sys- 
tematic treatment in an address is impossible. I shall, therefore, 
restrict my report to some phases of the problem to which I 
have made contributions, and to such papers as deal with re- 
lated questions and which have come to my notice. 

I sincerely appreciate this opportunity to bring together the 
results of several articles, all dealing more or less directly with 
the problem of complex roots, but scattered in various journals 
over a period of many years. 

Our equations are always polynomials, equated to zero, and, 
unless otherwise stated, with real coefficients. The notation 
fn(z) =0 indicates that the equation is of degree m. Merely to an- 
swer the question whether equations of high degree are ever 
actually solved, it may be pointed out that a table of natural 
sines and cosines, say for every 1’’, is nothing but a complete 
tabulation of the solutions of the equation x1!2%°— 1 =0, 


2. On Equations with Roots e®. We have in the literature many 
theorems on equations with roots of absolute value unity. One 
of the most interesting of these is Kronecker’s theorem: 

If the coefficients of an equation are ordinary integers, if the 
coefficient of the highest power is unity, and if all roots are of ab- 
solute value unity, then all the roots are roots of unity, and the 
equation 1s therefore solvable by radicals. 

A companion theorem is as follows. 

With the same restrictions on the coefficients, tf all roots are real 
and of absolute value <2, the roots are all of the form 2 cos rk, 
where k is rational. 

In two short notest I have derived necessary and sufficient 
conditions that an equation have some, or all, of its roots of 
the form e*. 


* An address presented to the Society, by invitation of the program com- 
mittee, Lincoln, Nebraska, November 30, 1934. 

{¢ Kempner, Archiv der Mathematik und Physik, (3), vol. 25 (1916), pp. 
236-242. Kempner, Téhoku Mathematical Journal, vol. 10 (1916), pp. 115- 
417. 
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THEOREM 1. To test any equation f,,(z) =0 with real coefficients 
for roots of the form e*, form the following polynomial, which has 
real coefficients and which contains only even powers of z: 


g2n(Z) (z + 1) f ) ¥,(u). 


The equation f,(z)=0 has among its roots some of the form e® 
when and only when p,(u) =0 has some positive roots. A necessary 
and sufficient condition ihat all roots of f,(z)=0 be of the form 
e” is that ,(u) =0 have only positive roots. 


For proof we have only to remember that z= (z;+7) /(z,—7) 
transforms the unit circle in the z plane into the axis of reals 
in the z, plane. The theorem may be extended in an obvious 
fashion to equations with complex coefficients, by considering 
f(z) -f(z) = F(z) in place of f(z). It is clear that the theorem ap- 
plies to much wider classes of functions than polynomials. If 
v,(u) =0 has only real positive roots, and if, in addition, the 
coefficients of f,(z)=0 are real integers and the coefficient of 
the highest power is unity, the roots of f,(z)=0 are roots of 
unity; and, by Kronecker’s theorem, the equation f,(z) =0 is 
solvable by radicals. 

Making use of the inverse transformation 2, =7(z+1)/(z—1) 
and retracing our steps, we now start from a given equation 
on(Z) =ao2"+a12""!+ - - - +a,=0 whose roots are all real, and 
form the new polynomial 


1 1 


z—1 z—1 


If @,(z) has integral coefficients, the coefficients of W2,(z) are 
also integers. It is easily seen that the coefficient of 2?" in 
say. If y=1, Yen(z) =0 satisfies the conditions of Kronecker’s 
theorem. Hence yY2,(z) =0 is solvable by radicals, and hence also 
¢2(z) =0 is solvable by radicals. 


THEOREM 2. Suppose that --- +a,=0 
has real integral coefficients and that all its roots are real. If 
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(@o—d2+a,—+ ---)?=1, the equa- 
tion is solvable by radicals.* 


3. On Equations with Roots a+pe®, where a and p are Rational. 
Of greater interest is the investigation of irreducible equations 
admitting roots of the form a+ pe, where a and p are both ra- 
tional.j Hand in hand with these equations go those admitting 
roots of the form a+ pe**, where a and p? are rational. 

We first characterize geometrically the set of points in the 
plane of complex numbers corresponding to the set of numbers 
a+pe”, where a and p are rational. The numbers are repre- 
sented by the doubly infinite set of circles whose centers are at 
rational points on the axis of reals, and whose radii are rational. 
We shall call this set of circles, or the corresponding set of num- 
bers, the set S. 

Similarly, we shall represent the set of numbers a+ pe*, where 
a and p? are rational, by a corresponding doubly infinite set of 
circles, and we shall call this set of circles, or the set of numbers, 
the set C. 

We shall consider the class of irreducible equations that have 
any complex roots on any circle S. For simplicity, unless other- 
wise stated, we shall assume that the coefficients are integral 
and that the equation is irreducible in the natural domain. This 
class consists of: 

(1) all irreducible equations f(z) =0 that have any (not neces- 
sarily all) complex roots of absolute value 1; 

(2) all equations derived from (1) by subjecting f(z) to a non- 
singular linear transformation z= (az,+))/(cz:+d) whose coeffi- 
cients are real and rational. 

In particular, all irreducible quadratic equations with ra- 
tional coefficients and negative discriminant belong to this class. 
The case (1) contains as a subclass all cyclotomic equations of 
prime degree. 

It turns out that if an irreducible equation with rational co- 
efficients has any root on a circle S, or C, the roots are sym- 


* A systematic examination of applications of Kronecker’s theorem may 
seem desirable. Compare also D. H. Lehmer’s article referred to at the end of 
the next section. 

t Kempner, On irreducible equations admitting roots of the form atpe®,---, 
Tohoku Mathematical Journal, vol. 13 (1918), pp. 253-265. 
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metrically distributed with respect to this circle. In this state- 
ment, two points are called symmetric with respect to a circle 
when each is the image of the other under transformation by 
reciprocal radii. This represents a natural extension of the sym- 
metry property of the axis of reals for all equations with real 
coefficients. 


THEOREM 3. When an irreducible equation has a complex root 
of the form a+pe, where a and p* are rational and where p is 
positive, but may be either rational or irrational, the roots of the 
equation are distributed in the following manner: 

(1) Besides a+pe**® there may be other pairs of complex roots 
with the same a and the same p, that ts, of the form a+pe*. 

(2) If cis a real root, p?/a is a root. 

(3) Any complex root not contained in (1) may, of course, be 
written in the form a+re**, with the same a. Then the four num- 
bers a+ret**, a+(p?/r)e*** are all contained among the roots of 
the equation. 


The complete extension of the symmetry property of the axis 
of reals, insofar as I can see, may be stated as follows. Let 
f(s) =2"+a,2""'+ --- +a,=0 be an equation with real, but 
otherwise arbitrary, coefficients, irreducible in the domain 
R(a,, - and let us assume that f(z) =0 has a root 
of the form a+pe”, where a and p? are both rationally expres- 
sible in terms of a;, - - - , @, with rational numerical coefficients. 
Then all roots of f(z) =0 are distributed in the plane of complex 
numbers symmetrically to the axis of reals and to the circle 
a+pe”. We may state also the following theorems. 


THEOREM 4. When an irreducible equation (in R(1)) with real 
coefficients has any complex root on any circle C, or on any straight 
line parallel to the axis of imaginaries and at a rational distance 
from it, the roots are symmetrically distributed with respect to this 
circle or straight line. 


THEOREM 5. The points of intersection of two circles S, which 
form a set that is everywhere dense in the complex plane, are exactly 
the set of roots of the class of irreducible equations az*+bz+c=0, 
where a, b, and c are rational and b? —4ac <0. In other words, the 
set of points of intersection of S are exactly the set of all complex 


— 
—- 
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irrationalities. The real quadratic irrationalities are given by the 
points of the axis of reals belonging to C, but not to S. 


It is easily shown that the complex (not real) points of the 
set of intersections of circles of the set C are the same as the 
complex (not real) points of the set of intersections of the set S. 
Through each point of intersection of two circles S (or C) pass 
an infinite number of circles S (or C). 

Therefore, any imaginary quadratic irrationality may be 
represented in an infinite number of ways by the form a+pe*, 
where a and p are rational, but any root of an irreducible equa- 
tion of degree »>2 can be represented in at most one way 
in this form. Consequently, for >2, the equation ai+p.e™ 
=a2+pee™, where aj, pi, and pe are rational and where p; 
and pe are positive, implies that a; =a2 and that p: =pe. 

It follows from the preceding theorems that if an irreducible 
equation of degree n >2 has roots on a circle S, these roots can- 
not be among the points of intersection of circles S. 


THEOREM 6. An irreducible equation cannot have roots on more 
than one circle S unless the circles intersect. On every circle S lie 
roots of irreducible equations of degree n>2, some of whose other 
roots lie on tntersecting circles. 


THEOREM 7. An irreducible equation f(z) =0 can have two roots 
ayt+pie™, where ay, a2, pi, and are all rational and 
a2, only when | p1—pe| <pit/pe. 


THEOREM 8. When an irreducible equation f(z) =0 has complex 
roots of the form z=a+pe”, where a and p are rational, a and p 
may be found by a finite number of rational operations and the 
extraction of a square root. Therefore, if such an equation has a 
root on any circle S, this circle, but not necessarily the root itself, 
can be constructed by ruler and compasses. 


Let us give a numerical illustration. The irreducible equation 
1624+ 722*+ 962? + 602+ 15 =0 is found to have roots of the form 
a+pe”, wherea = —1/2,p=1/2. The equation has two real roots, 
(—(21)!/2—9+ and on the circle —1/2 
+e%/2, two complex roots, { (21)!/2—9+ [30— 10(21)1/2]1/2} /8. 

Necessary and sufficient conditions for an irreducible equa- 
tion to have complex roots of the form a+pe”, where a and p 


= 
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are both rational, are obtained by extending our theorem on 
necessary and sufficient conditions that an equation have roots 
of the form e*. 


THEOREM 9. An irreducible equation f,(z)=0 has among its 
complex roots some roots of the form a+pe”, where a and p are 
rational, when and only when both of the following conditions are 
satisfied. 

1. The system of equations 


= f(a), *(maga + a1) = f’(a) 


has a rational solution (a, p). 
2. Letting z=a+pz and forming the expressions 


g(z1) = p21), 

= (6:2 + e( = ¥a(st) = 
2-1 


we see that the equation y,(u) =0 must admit at least one positive 
root. 


As applications, we shall mention the following cases. 

(1) Klein’s famous Modul-figur, consists, in the plane of com- 
plex numbers, of the set of circles of radius 1/2 with centers on 
the axis of reals at 0, +1, +2, - - - , and the system of straight 
lines parallel to the axes of imaginaries which are tangent to 
these circles and the system obtained by continued reflection 
of this figure with respect to itself by reciprocal radii. We may 
state the following property. An irreducible equation cannot have 
roots on more than one circle of the Klein Modul-figur. 

(2) Every irreducible equation of degree <10 with a root on any 
circle S is completely solvable by radicals.* 


(3) When an irreducible equation aoz"+a,2""!+ -- - +a,=0 
has a complex root of the form a+pe, where a and p are rational, 
the equation ayy" --- +a, has a rational solu- 
tion (x, y). 

When a and p are both integers, the diophantine equation 
doy" =aox"+ -- + +a, has an integral solution (x, y). 

If the equation aoz"+a;2""!+ - - - +a,=0 has a complex root 


* This can probably be extended to give a result corresponding to the last 
theorem of §2. 
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on a circle of rational radius touching the axis of imaginaries at the 
origin, the equation a,x""!+<a2x"-*+ --- +a,=0 has a rational 
root. The equation 162*+722°+96z?+60z+15=0, which was 
mentioned above, comes under this heading. The equation 
72x? +96x?+60x+15=0 has a root —1/2. 

In connection with this work I may mention a paper by D. H. 
Lehmer.* He defines a polynomial as quasi-cyclotomic if all zeros 
lie on a circle of radius r, with center at the origin, and if the 
arguments of the roots are commensurable with 27. His work 
is connected with our own by his Theorem 3, which is as follows. 

Let the roots of f=0 lie on the circle of radius r. Then r is an 
integer or the square root of an integer, according as n is odd or even. 

An interesting result is his Theorem 6, which reads: Let the 
n=2k roots of f(x) =0 be of the form (R)'?-«€,, where € are roots 
of unity, and where R has no square factor. Then 


(R)'?2 (2k)! 
(k + 1)\(k — 1)! 


4. On the Distribution of Complex Algebraic Integers of a Given 
Degree. The identification of the points of intersection of the set 
of circles S (or C) with the totality of complex quadratic irra- 
tionalities (integral or otherwise) suggests the problem of in- 
vestigating the distribution in the plane of complex numbers of 
the totality of algebraic integers x+1y of a given degree n=2, or 
n=3, and so on. 

The problem appears artificial on account of the fact that the 
totalities considered do not form domains of rationality. Yet, 
it seems to possess a degree of interest. For m =2, the solution is 
given practically by inspection. In the (x+7y) plane, these 
algebraic integers are obviously distributed along the lines 


x=k/2, (k=0, +1, +2,---). The distribution is the same 
for all lines x=0, +1, +2, ---, and is also the same for all 
lines x= +1/2, +3/2,---, but it differs in the two sets. The 


density increases with increasing |y|, and so on. In addition, 
the integral real quadratic irrationalities are of course distrib- 
uted everywhere densely on the axis of reals. 

For the set of all cubic algebraic integers (and eo ipso those 


* Quasi-cyclotomic polynomials, American Mathematical Monthly, vol. 39 
(1932), pp. 383-389. 
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of higher order), the first suspicion is that they form in the plane 
of complex numbers an everywhere dense set. Even if this 
should be so, there is some indication that an examination of 
the fine structure of the distribution with respect to “gleich- 
miassige Dichtigkeit” may be of interest. 


5. On the Use of the Arguments of the Coefficients for Separating 
Complex Roots. In a previous paper,* I emphasized the possi- 
bility of using systematically the arguments of the coefficients 
of an algebraic equation for the separation of the complex roots. 
Descartes’ rule for an upper bound of positive (negative) roots 
is the first, very special, case of theorems of this kind. I used a 
very obvious, yet rather powerful, tool. Consider the equation 


f(2) = a,2"+---+a=0, 
a, = r,e%, 0 < 0, < 2x, r, > 0; = pe**, < 2x, p> 0. 


In the plane of complex numbers, mark from the origin the 
half-rays or vectors 0,. The vector 4 is to remain in its original 
position. The vector 4; is to rotate in a positive sense with a 
constant angular velocity w, while for 1=1, 2, - - - , m, the vec- 
tor 6, rotates with a uniform angular velocity yu times that of 4. 
Vectors 0, for which r,=0, that is, for which the coefficient 
a, =0, are to be ignored. At any moment the vectors give the di- 
rections of the vectors representing the terms a,2" = 17, p*e*@#t#9) 
These directions depend only on ¢ and the 6, of the coefficients. 
Most of the theorems are immediate consequences of the fact 
that the sum of vectors from a common point can certainly not 
vanish when it is possible to draw a line through the point such 
that all vectors lie on one side of the line. This fact is used, for 
example, in the simplest proof of the Gauss-Lucas theorem that 
all roots of f’(z) =0 lie inside or on the boundary of the smallest 
convex polygon around the points representing the roots of 
f(z) =0. 

It is, therefore, impossible to have roots of f(z)=0 in any 
sector (vertex at the origin) a<@<§ for which all vectors lie 
on one side of a straight line through the origin. Such ¢-inter- 
vals—if they exist—are determined by simple inequalities, 


* Kempner, Uber die Separation komplexer Wurzeln, Hilbert Festschrift, 
1922, pp. 49-59 (= Mathematische Annalen, vol. 85 (1922), pp. 49-59). 
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which in many cases enable us to determine sectors, depending 
only on the arguments of the coefficients, not on the absolute 
values, and which are free of roots. 

The efficiency of the method is a function of the number of 
coefficients different from 0, rather than of the degree of the 
equation. For trinomial equations, a complete separation of the 
roots is given; for quadrinomial equations, and m-nomial equa- 
tions, m>4, valuable information usually is obtained, even 
though complete separation may not be accomplished. 

The results are indicated by the following illustrations. 

(1) Binomial equations: a,2"—ay=0. We obtain n null sec- 
tors (half-rays), each of which contains exactly one root, in 
agreement with elementary algebra. 

(2) Cubic equations: 2*+az+b=0, where a and b are real. The 
three roots have arguments 41, ¢2, 63 (6 =0, real positive roots; 
o=m7, real negative roots); 


5x 


5x 


(3) Trinomial equations: 
Om ao=roe™. For such equations, we obtain com- 
plete separation of roots, one in each of m sectors.* Special 
precautions are required for dealing with the cases in which 
two sectors adjoin and the roots lie on the common boundary 
line. As an illustration, consider equations of the form 


a5z5 + aoz? + ap = 0, 6, = 0, 6. = = 57/6. 
Each of the following five sectors contains exactly one root: 
66° < 84°, 129° < 138°, 204° <3 < 210°, 282° 309°, 


* Partly overlapping results by Nekrassoff, Mathematische Annalen, vol. 
29 (1887), pp. 413-440, who has different aims and uses entirely different meth- 
ods. 
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324° <5 < 354°. The openings of the sectors are, respectively, 
18°, 9°, 6°, 27°, 30°. This result holds, of course, independently 
of the absolute values of the coefficients. 

(4) Quadrinomial equations. As an illustration, consider any 
equation of the form 


+ + a525 + ao = O, (do, @s, Ay, > O). 


Roots are separated into ten sectors: 18°<¢,:<36°, 36°<¢» 
<60°, 90°<¢;<108°, 108°<$,<140°, 162°<¢;<180, 180° 
< 198°, - - - , 324° << 342°. The largest opening is 32°, 
the smallest 18°. Complete separation is partly due to the high 
degree of the equation, despite the small number of terms. For 
quadrinomial equations, complete separation is no longer guar- 
anteed. 

(5) Mechanical Determination of Sectors. An instrument of the 
type of a “planetarium” which permits m vectors to rotate from 
arbitrary initial positions with angular velocities which are to 
each other in the ratios m:n—1: - - - :3:2:1 would permit the 
mechanical determination of sectors free from roots. An ordi- 
nary watch will separate all roots of any equation of the form 


22+ (a>0,6>0,A =1ori11). 


Consider, for example, any equation of the form 2” +a,z+a)=0, 
(do, >0). We obtain 
° ° ° 


4 1 9 
15°<gi< 1675 4955 » 8155 2 


° 


6° 3 
105° 11455 , 135°<gi< 1475 , 180°, 


with a similar distribution of 6 sectors from 180°—-360°. The 
openings of the sectors are, respectively, 
N. Obreschkoff* has given another extension of Descartes’ 
rule. His result is as follows. 


* Sur un probléme de Laguerre, Comptes Rendus, vol. 177 (1923), pp. 102- 
104. 
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Let f(z) =ao+aiz+ -- - +a,2"=0 be an equation with real co- 
efficients. The number v, of variations in the sequence of its coeffi- 
cients is not only an upper bound for the number N of its positive 
roots, but also an upper bound for the number of roots lying in the 
sector larg 2| <a/n. If the number of roots lying in the sector is 
less than vq, 1t 1s smaller by an even number. 

I. J. Schoenberg was kind enough to let me see the manu- 
script of an article soon to be published in the Annals of Mathe- 
matics. In the first part of this paper, the author introduces an 
ingenious extension of the notion of a sequence of variations 
in signs to the case of complex coefficients. This is accomplished 
by the so-called separating double-sectors for the coefficient ap. 
To define these, we imagine in the plane of complex numbers 
the coefficients a@, marked by their corresponding points. 
Through the origin draw any two straight lines J; and /, that 
divide the plane into four consecutive sectors Ai, Az, Az, As, 
with the following property: one of two pairs of opposite sec- 
tors, for example, A; and A3, shall not contain any of the points 
a, in their interior. If the common opening of A; and A; is 
v,(0<WS7),we call A; and A; together a separating double-sector 
(of aperture y) for the coefficients a,. In going in order through 
the sequence of points do, a1, - - - ,@,, we count the number of 
times we have to pass from the sector Az to A, or vice-versa. 
This number is called v,(.S) and is equal to the number of varia- 
tions of the sequence do, - - - , @, with respect to the separating 
double-sector S=(A,, A3). It is evident that for all a, real, v.(S) 
(taking Y=) agrees with the ordinary number of variations. 
Schoenberg then derives the following theorems. 

Let f(z) =ao+aiz+ -- - +a,2"=0, where a, is real or complex. 
In the plane of complex numbers, mark the points a, and draw a 
separating double-sector of aperture py, (0<WS7), as can always be 
done. Then the number of roots lying in the sector larg 2| <y/n is 
either equal to va(S), or ts less by an even number. To find an 
upper bound for the number of roots of f(z) =0 lying in a sector 
whose bisectrix makes with the positive axis the angle 0, a substitu- 
tion 2=e'z' is carried out. 

For any sector with the vertex at the origin and of aperture 
2x/[n(n+1)] an upper bound of the number of roots lying in the 
sector can be established. 

The rest of this interesting paper deals with applications of 
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these ideas to equations of various types which fall outside the 
frame of this paper. 


6. Absolute Value, Argument, Real Part, Pure Imaginary 
Part of Complex Roots. 


Given 
f(z) = + + + a, = 0, 


where z=a+78 =re’® is a root. It is well known that, if a, are 
integers, then 7, cos ¢, sin ¢, a, and 8 all are algebraic numbers, 
and satisfy equations with integral coefficients which may be de- 
rived from f(z) =0 by rational operations. It would be of some 
interest to have a systematic examination of these equations 
and their use for the determination of complex roots. At present 
the information is badly scattered in the literature.* 

Difficulties arise in all cases: (1) because the equation for 
any of r, cos ¢, sin ¢, a, B is of high degree, usually n(m—1), 
with reductions in some cases for real rational a,: for example, 
for aoz?+ +a2=0, z=a+i8, B satisfies a second degree 
equation, but a satisfies a linear equation, with integral coeffi- 
cients; (2) because the equations necessarily have a large num- 
ber of extraneous roots, which are sometimes not readily distin- 
guishable from the roots required. Some aspects of a systematic 
examination and classification of methods of solution along 
these lines are contained in two master’s theses. f 

It is interesting to note that the equations for r and for e** 
are related in structure (for a, real). Let z =re® =rs be a complex 
root of f(z) =0. 

Then r/s is also a root, and hence 


+ ayr™ 14+ --- + a,17m5 +a, 


0, 
0. 


Eliminating, we have as equations for 7 and s, respectively, 


* References are too numerous to mention. We refer to just one interesting 
paper by Hayashi, Téhoku Mathematical Journal, vol. 3 (1913), pp. 110-115, 
for the derivation of the equation for the absolute values. The degree of the 
resulting equation is incorrectly given [n, instead of n(n —1)]. 

7 R. C. Huffer, University of Illinois, 1920, and A. W. Randall, University 
of Colorado, 1929. 
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In connection with the problem of numerical solution of 
algebraic equations, I should like to mention a doctoral thesis* 
by A. J. Lewis on the solution of algebraic equations by infinite 
series, since it has not yet been published. As in practically all 
theoretical work on solutions by series, the background is given 
by Lagrange’s expansion. The method was elaborated by 
Nekrassoff (1887), McClintock (1894), Lambert (1908), and 
others. However, the whole work is very fragmentary in char- 
acter; and Lambert’s work, in particular, is open to criticism. 
Lewis has given a more detailed treatment of the expansions 
than was contained in the literature, has corrected Lambert’s er- 
rors, and has made a step in advance in the investigation of the 
convergence conditions. Trinomial and quadrinomial equations 
are discussed in detail, and the treatment for equations with 
more than four terms is indicated. The formulas are very com- 
plex, but this is due to the intrinsic difficulties of the problem. 


7. On a Graphical Method for the Separation and Computation 
of Complex Roots. The next article I shall take the liberty of 
reporting more in detail, since it is published in a journal not 
easily available to mathematicians.f In this paper, three meth- 
ods of graphical representation are compared.{ 

(1) The representation familiar from conformal mapping, or 
from proofs of the Fundamental Theorem of Algebra: f(z) 
=f(x+ty) =u(x, y)+10(x, y). The intersections of the two 
curves u(x,y)=0, v(x,y)=0 in the xy plane give the roots 
x+zy of f(z) =0. In order to find the number of roots of absolute 
value p, where p is arbitrarily assigned, it is only necessary to 


* A. J. Lewis, Thesis, University of Colorado, 1932. 

{t Kempner, On the separation and computation of complex roots of algebraic 
equations, The University of Colorado Studies, vol. 16 (1928), pp. 75-87. 

t For the careful computation of all graphs of this paper, I am under great 
obligation to my colleague Professor C. A. Hutchinson. 
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find the number of points of intersection of the u curve and the 
v curve inside a circle about the origin of radius p. 

It is possible, by means of Kronecker’s characteristics, or by 
equivalent methods, to determine from the figure the number of 
roots of f(z) =0 inside any closed contour by considering nothing 
but the points on the contour where the curves u=0 and v=0 
enter and leave the region. We shall use this idea later. 

(2) Let 

S(@) = + --- +o =0, 


where (cos yeti sin yx) =rpe*%, and z=x+1ty 
=p (cos ¢+1sin d) =pe*®. 


We plot 
“u= r.p*cos (ko + r.p*sin (kb + vx) 
k=1 k=1 


in the following manner. For p, the absolute value of z, we select 
an assigned value which is held fixed. Then u and v become func- 
tions of ¢ alone, say u=y(¢), v=y2(o). Thus u and v may be 
considered as coordinates of a curve given in parametric form, 
with ¢ as the parameter. We plot the curve in a uv system of 
coordinates, choosing the horizontal axis as the u axis and the 
vertical axis as the v axis. 

(3) We now plot the two curves u =y(¢), v=y2(¢) separately, 
but both in the same system of coordinates. For both curves, 
we take the horizontal axis as the ¢ axis, the vertical axis as the 
u axis for the first curve, and as v axis for the second curve. 

Discussion of (2). To a circle in the complex z plane about the 
origin, of given radius p, corresponds in the uv plane, which we 
may also interpret as the (w+iv) plane, a certain closed curve. If 
N is the number of roots of ¥(z) =0 in the circle, and if we assume 
that no roots lie on the circumference, each root being counted 
with its proper multiplicity, then 27N represents, as is known 
from the theory of analytic functions of a complex variable, the 
change of amplitude when the curve is once completely traced 
in the (u+7v) plane. 

Simple considerations of a type familiar from the theory of 
characteristics will permit us to read off this number N from our 
representation (2). To find the number of roots in a circle of 
radius p about the origin in the z plane, we consider the cor- 
responding closed curve (2). Trace the curve (in either the posi- 


é 
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tive or negative sense), starting from an arbitrarily chosen point, 
and mark as points 6 all points where the curve crosses the u 
axis (v=0), and as points a all points where it crosses the v 
axis (~=0). Now we determine the change in the argument be- 
tween two consecutive points of intersection with the wu axis, 
as we move along the curve. We call this argument y, since it is 
not necessarily the same as the parameter ¢, the argument of 
the construction (3). The simplest case occurs when our points 
a and 8 form an alternating sequence, aBaB--- a8. To a 
change from any a to the next a corresponds each time a change 
of z in y. All these changes are in the same sense; that is, y is 
either each time increased by 7 or each time decreased by 7. 
Hence y increases or decreases by mz, where m is the number of 
B, and there are therefore m/2 roots in the interior of the circle 
of radius p about the origin in the (x +7y) plane. But the number 
m=2N is correctly expressed by the following theorem for the 


special case --- =q=1,v=m. 
THEOREM 10. Denote by B:1, Be, - - - the points of intersection of 
curve (2) with the u axis, by a, Q2,--- the points of intersection 


with the v axis, and consider the chain of the points a and B cycli- 
cally closed. If then, starting from any a, the number of B points 
between two consecutive a points is, respectively, &, €2,°-°-*, Ek; 
then the number of roots z of f(z) =0 for which | s| <p ts given by 
the formula 


2N = 1) (ect) + let) 


| k=l 


I omit the proof that this formula, derived for the special 
case of alternating a and 8, is true for any sequence of the 
a and B. These considerations lead to the following rule. 


RULE. Write down, in proper order, the numbers which repre- 
sent the number of B in the successive a-intervals. Under the first of 
these numbers write +1, if it is odd, or —1, if it is even. Write 
then under each odd number that one of the two numbers, +1 or 
—1, which is already written under the last preceding number; but 
under each even number write +1 when —1 is written under the 
last preceding number, and write —1 when +1 was written under 
the preceding number. The sum of all of these numbers +1, —1 is 
necessarily even, and the absolute value of the sum is 2N. 


| 
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Discussion of (3). From the constructions (2) and (3), we 
find the following relation between them. 

It is clear that to every curve (3) will correspond uniquely a 
curve (2), since for every value of ¢ both u and v are uniquely 
determined. On the other hand, from a curve (2) a curve (3) is 
not uniquely determined : according to the manner of parametric 
representation of u and v by means of ¢, an infinite number of 
curves (3) may be obtained. All curves (3), however, which cor- 
respond to the same curve (2), are related to one another in a 
simple fashion. 

The theorem is applied, without modification, to (3). In our 
curves (1), drawn for a given value p, and with a common @ 
axis, the intersections of the u curve with the ¢ axis and the 
intersections of the v curve with the @ axis correspond in a 
definite manner to the intersections of the uv curve in (2) with 
the u axis and the v axis. To make this clear, it is only necessary 
to consider carefully a special case. Consider the uv curve 
treated in Fig. 12.* To this figure corresponds (schematically) 
Fig. 2; as is immediately checked. In comparing these two 
figures we have only to keep in mind that every time the curve 
(2) crosses the u axis, we have v=0, so that in (3) the v curve 
crosses the ¢ axis, and that every time the curve (2) crosses the 
v axis, we have u=0, so that in (3), the u curve crosses the 
@ axis, and finally, that the combined sequence of the points 
of intersection of the u curve and the v curve with the ¢ axis 
follow in the same order as the af sequence in (2). 

It follows that we get exactly the same a6 sequence, whether 
we read it off from (2) or from the sequence obtained by marking 
in (3) as @ all points of intersection of the u curve with the ¢ 
axis, and as 8 all points of intersection of the v curve with the 
@ axis. We have, therefore, the following theorem. 


THEOREM 11. Plot, as explained under (3), the curves 
u = cos (kp + v= sin (ko + 14), 
k=0 k=0 


for a given p. Name the points of intersection of the u-curve with 
the axis, 04, Q2,--- , and name the points of intersection of the 
v curve with the axis, Bo, Consider the combined se- 


* The index 2 in 12 indicates that the figure refers to construction (2), etc. 
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quence of the a and B points in their natural order of magnitude. 
Assume the sequence closed cyclically, and let &, €,--- , €, re- 
spectively, be the number of B points between two consecutive a 
points; then the number N of roots of f(z) =0 for which | 2| <p is 
given by the formula 


2N =| (— 1) | 
k=1 


This number may be evaluated by the rule stated under Theorem 10. 


EXAMPLE. The following example is only schematic, and 
shows how the rule and the theorem are applied to more com- 
plicated cases. 

Figures 1; and 2; are corresponding figures. From either of 
them we read off the a8 sequence, and apply the rule or the 
theorem. 


Baa B aaa BBa Ba Ba BB 
0 1 0 0 Z 1 1 3 
—1 —1 +1-1 +1 +1 +1 +1 
Using the rule, we find 2N =2, N=1. 
Using the theorem, we find 
e+1= 1 2 ee 3 2 2 4 
+(—1)s=2, N=1. 


If we had started with 6 intervals and had considered the num- 
ber of a in each B interval, we should have found 


Baa Baaa BB a Bp a Ba B B 
€.= 2 3 0 1 1 1 0 0 
e+1= 3 4 1 2 2 2 : ae | 


+(-1)*=2, W=1. 

A pplication of Theorem 11. We make use of the theorem to 

isolate the complex roots of an equation in concentric rings 


(p1<p<p2) around the origin as center. The procedure is as 
follows. For two given positive values of p, p2>p1, we assume 
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that the curves u(¢) and v(¢) are drawn in our system (3). By 
Theorem 10 we read off the numbers N(,), N(p2) of roots of 
absolute value smaller than p; and pe, respectively. Then N(p2) 
—N(p:) gives the exact number of roots in the ring-region 
p2>p>p.. Here we assume, since this case is easily disposed of, 
that no root lies on the circumference of either of the circles 
OF pe. 

In case we have the same sequence of a and 6 for p; as we have 
for p2, there are obviously no roots in the ring-shaped region 
| </pe. 

Since there exist well known simple upper and lower bounds 
for the absolute value of the roots of a given equation, we may 
start out with two such values for p, and, by appropriate choice 
of new values of p, we can separate our roots according to their 
absolute value. 

It is a simple matter to plot curves of the kind used here, 
particularly because we are really interested only in the se- 


Fic. 3; 
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quence of the points of intersection of the u curve and the v 
curve with the @ axis. The chief advantage of the present 
method may lie in the fact that instruments exist with which 
curves of type (3) can be drawn. Such an instrument is the 
Michelson-Stratton harmonic analyzer, for which an accuracy 
of about 1 per cent to 2 per cent is claimed. 

EXAMPLE. Consider the equation 


227+ iz+1=0,u =p?cos2¢ — psing + 1,v =p’? sin 2¢ + pcos 


Figure 3, gives the classical construction, which may be used for 
the verification of our results. Figure 42 gives, for p=1/2, the 
parametric representation of u=u(p), v=v(@) in a uv-system of 
coordinates, while Fig. 53 gives u and v both as functions of ¢. 
Figures 6. and 7; are the corresponding curves for p =1; Figs. 
8 and 9; are the curves for p =3/2; and finally Figs. 102 and 11; 
are the curves for p=2. 
For p=1/2, from Fig. 53, we have 


BB 
«= 00, 2N=(—1)!+(-2)?=0, N=0. 


This is obvious, since the u curve does not intersect the @ axis 
at all. For p=1, from 73, we find 


aBaBBB 


Il 


For p=3/2, from 93, we get 
aBaBBB, N =1, asforp =1. 
For p=2, from 113, we obtain 
aBaBaBpaBp 
e= 1 1 1 1 


Therefore, the equation z?+iz+ 1 =0 has one root in the ring- 
region 1/2<p<1 and one in the region 3/2 <p <2. This agrees 
with Fig. 3;. The absolute values of the two roots are 5/2 
—1/2=0.6+ and 5'?/2+1/2=1.6+, respectively. 

It will be recognized immediately that the graphs shown in 


— 
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Figs. 53, 73, 93, 113; also give information concerning the argu- 
ments of the complex roots. 


8. Separation of Roots in Sectors. Without insisting at the 
present time, I shall only indicate how these same ideas can be 
utilized to separate roots in sectors from the origin by inter- 
changing the role of r and ¢ in the former process. 


FIG. 62 


Given the equation f(z) =aoz"+ - - - +a,=0, where @, is real 
or complex, and again z=pe**=p(cos ¢+i sin @). Now, for a 
given ¢, we have u=u(p),v=v(p). Plot u(p) and v(p) in the same 
system of coordinates, with p as horizontal axis, and with u 
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and v as vertical axes. We know how to find upper and lower 
bounds p; and p» for the absolute value of the roots. The value of 
p (always >0) is to be restricted to the interval between these 
upper and lower bounds. If the sequence of a and 8 is the same 
in the interval p2) for =¢; and for there are no roots 
of f(z) =0 in the sector ¢iS ¢Z¢@2. If the sequence of a and 8 
is not the same for 6=¢; and @=q@», there are roots in the 
sector and the exact number may be read off by the same kind 


Fic. 73 


of theorem or rule as that which was derived above. The details 
of this are reserved for another occasion. At present, we restrict 
ourselves to the verification that 23—2z—2=0 has a root in 
the sector 37/4 <¢ Xz, with 1<p<3/2: 

2—2z—2=0. 


u = p*? cos 36 — 2p 2, v= p* sin — 2psin ¢. 


1 1 
= 32/4: (2)1/2p8 (2)!/2p —2, v (2)"/2p?— (2)!/2p; 


— p?+ 2p — 2, v= 


6 
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Since for ¢=7, v reduces to 0 identically, we select ¢=7—e, 
where € denotes some small positive angle. Then, from the rela- 
tion lim,.o[(sin €)/e]=1, we find v=3p%e—2pe, approximately. 

For ¢=3r/4, our u curve, for 1<p<3/2, lies entirely above 
the p axis, and the v curve crosses the p axis once (no a, one £). 


-3 18, -2 —1 0| \& 1 2 3 Ba\Be 


Fic. 82 


For ¢=z7—e, neither the u curve nor the v curve cut the p 
axis for 1<p<3/2 (noa, no 8). There must be roots of f(z) =0 
that lie in the sector 37/4 <@ <r. 

We do not have, as we did before, an instrument available for 
drawing the curves, but, on the other hand, they are of par- 


v 
-2 
\ 
—3 
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ticularly simple type, namely, polynomial curves. Besides we 
need, of course, not the curves themselves, but only their (ap- 
proximate) intersections with the p axis. 


9. Polynomial Curves and Complex Roots. I shall discuss one 
last topic also with some degree of completeness, in order to 
piece together ideas which are scattered and more frequently im- 
plied than explicitly announced in the literature. 


Fic. 9; 


The importance of polynomial curves for the real roots of 
algebraic equations is fundamental and obvious. But the sys- 
tematic use of graphs for the purpose of locating complex roots 
of equations of the form f(z) =ao2"+ - - - +a,=0, where a, 
is assumed, for simplicity only, to be real, is not so widely 
known as it should be, except for small values of , such as oc- 
cur in the problem of the location of complex roots of quadratic 
or cubic equations with real coefficients. 

The character of the real graph of f(z) is of great influence in 
some of the graphical methods for complex roots. I have clas- 
sified real polynomial curves of any given degree m according to 
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their geometric behavior, and have given existence proofs for 
the various types.* 

Of particular importance for our immediate purpose are cer- 
tain considerations which permit us to represent in a natural 


2 4 BAB 


Fic. 102 


manner some parts of the complex graph, and which, although 
the germ of the method goes back at least to Poinsot and his 


* Kempner, On the shape of polynomial curves, 1. Téhoku Mathematical 
Journal, vol. 37 (1933), pp. 347-362; Part II is to appear in the same journal. 
Compare also American Mathematical Monthly, vol. 40 (1933), pp. 469-470. 
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cordes idéaux, are not given the place they deserve in modern 
presentations.* 

The simple idea consists in an attempt to make the best of a 
bad situation. Since it is impossible to give in a three-dimen- 


3 
1 \ 
1 ! \ / 
6 : | ! \ / 
| \ 
\ / eer 
i / 
\ / 
/ 
\ / 
/ 
4 
\ 
Fic. 113 


sional space a simple satisfactory representation of a function 
w =f(z) of a complex variable, we try to save as much as we can. 
We may use our three dimensions as follows. Consider only real 
values of f(z), but permit z to assume all values, real and com- 
plex, which make f(z) real. 

Letting f(z) =f(x+iy) =u(x, y)+i0(x, y)=u+iv, we have, 


* Phillips and Beebe, Graphical Algebra, 2nd ed. (1887) (out of print), make 
much use of it in an elementary study of equations of the second, third, and 
fourth degrees with real coefficients, but with no attempt at a systematic treat- 
ment. In particular, the curves u(x, y)=0, v(x, y)=0 in the xy plane are not 
considered. The space curves in our text are to be found in Phillips and Beebe. 


/ \ 
\ \ 
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for f(z) real, v(x, y)=0, f(z) =u(x,y). If we use a horizontal 
xy plane for z=x-+17y, and a vertical axis for u(=f(z) real), 
we shall obtain a space curve in our (xyu)-system of coordinates. 
The projection of the space curve onto the xy plane consists 
exactly of the curve v(x, y) =0, and the functional value is the 
corresponding value of u(x, y). 

It is evident at once that this representation is related to the 
theory of functions of a complex variable, since the space curve 
is also given as the curve of contact between the surface of ab- 
solute values |f(z)| and the surface of real values u(x, y) for 
u(x, y)>0, and as the curve of contact between the surfaces 
~ | f(z) | and u(x, y) <0. Again, the space curves may be inter- 
preted as the intersections of the surface u(x, y) with the cylin- 
drical surface v(x, y) =0. 

These space curves serve as a skeleton for the function. They 
are important for us, (a) because we can study their general 
character with great ease, (b) because by means of them we can 
always represent all roots of the equation. 

(1) For a real function f(z) =ao2"+ --- +a, with maxima 
and minima of the first order (f’ =0, f’’ #0) we obtain a space 
curve of the following type. First we have in the xu plane the 
ordinary real graph, but in addition we have attached to each 
extreme a handle, symmetric to the plane of the real graph 
and penetrating it under a right angle. These handles always 
extend to infinity and their projections in the xy plane are (as 
stated) branches of the curve v(x, y) =0, the properties of which 
are well known, sometimes under the name of Lucas’ stelloids. 
In certain cases such a handle lies entirely in a plane, which 
then is for reasons of symmetry parallel to the yu plane. This 
happens, for example, for the real quadratic function* az?+ bz 
+c, in agreement with the fact that v(x, y) =0 now consists of 
a line parallel to the y axis, besides the x axis. Similarly, for real 
polynomial curves with a line of symmetry, which then neces- 
sarily have a maximum or a minimum on this line, and which 
we assume to be an extreme of the first order, the handle at- 
tached at that extreme lies entirely in the plane perpendicular 


* Leading to one of the classical constructions of complex roots of quadratic 
equations with real coefficients. 
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to the xu plane and contains the line of symmetry. For example,* 
consider the function f(z) = 24 — 22? = (z — (2)'/*) Here 
v(x, y) =0 consists of the curves y=0, x =0, x?—y?—1=0. The 
projection of the real part of the curve is y =0; the line x =0 con- 
tains the handle attached at the maximum; the branches of 
the hyperbola are the projections of the two handles attached 
at the minima. These last two handles therefore do not lie in 
planes. The distance of their points from the xy plane increases 
monotonically as the projection moves along the hyperbola to 
infinity. 

The same considerations may be applied to transcendental 
functions. For sin z, the result is particularly simple: v(x, y) =0 
consists in the xy plane of the axis of reals (the x axis) and of the 
equidistant lines x= +(2k+1)m/2, (R=0, 1, 2,---). Hence, 
all handles lie in planes perpendicular to the xu plane and 
parallel to the y axis. The handles are nothing but catenaries, 
attached at the maxima and minima, and extending upward 
(downward) from a maximum (minimum). Our later remarks 
concerning the utilization of these extended graphs for the 
location of complex roots of an equation ao2"+ --- +a,=0 
are immediately applicable to the location of complex roots of 
sin z=c, (c real, |c| >1). Applied to any two real circles in a 
plane, the handles lead to a neat interpretation of the circular 
points. 

For any equation of the form aoz"+ --- +a,=0 with real 
coefficients and with the maximum number n—1 of real ex- 
tremes (for example, f(z)=(g—ai;) --- (g—an)+B, 
~~ --- ~a, real, B real), the structure of the real curve and 
all of its handles is of such simple character that it seems a 
pity not to make use of it for approximating the location of 
complex roots. With but little practice, it is a simple matter to 
estimate roughly the points where this space curve intersects any 
plane parallel to the xy plane. These points, read off in the form 
x+ty, give a first approximation to all roots of f(z) +c=0. 


* The qualitative—not quantitative—aspects of many theorems in an in- 
teresting article by H. B. Mitchell, Transactions of this Society, vol. 19 (1918), 
pp. 43 ff., become intuitively evident by these geometric interpretations. 

See Fig. 17 for s*—22?+41. 
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I may mention that the Fundamental Theorem of Algebra 
can be given a very clear interpretation by means of such space 
curves, and that a proof of the theorem, essentially equivalent 
to the proof of the existence of m intersections of u(x, y) =0 and 
v(x, y) =0, can be constructed. 

(2) It is now necessary to free ourselves from the restriction 
that az"+ --- +a, have the maximum number of real ex- 
tremes. I can only indicate the state of affairs. 

To illustrate, we consider successively the functions* f(z) 
23—ez, 23, s3+e2, 0<e<1. As e—0, the two ex- 
tremes of z*—z move together, coalescing for «=0, and leading 
to a point of inflection with a horizontal tangent (f’(z) =0, 
f''(z) =0).f As we pass from z* to z*+2, the point of inflection is 
still existent, but the roots of f’(z) =0 have become imaginary. 
(For higher degree polynomials, the points of inflection gener- 
ated by the coalescence of consecutive maxima and minima 
may be ultimately lost.) To trace this in detail requires distinc- 
tion between natural points of inflection on polynomial curves, 
such as must exist between any two consecutive extremes, 
and accessory points of inflection, due to such coalescence of 
consecutive maxima and minima. At a natural point of inflec- 
tion, | f’(z) | is a maximum, at an accessory point, | f’(z) | is a mini- 
mum.){ Schematically, ~,— (maximum), (minimum). 

As € goes from 1 to 0 to —1, the handles first run together, 
and then completely separate from the real curve. The projec- 
tions of our space curves onto the xy plane are given by the cor- 
responding equations v(x, y) =0. 

Even if the accessory points of inflection should become lost 
(compare, for example, the succession of curves: z*—22?+1, 
z*—ez?+1, 24+1, 2t+es?+1, 24+22?+1) these branches are 
still present in the space curve. They can in no way be destroyed 
provided a)#0, but can only be changed in shape and position 
and in the manner in which they are connected with each other, 
however the coefficients a, may vary. 


* See Figs. 12-16 for e=1, 0.1, 0, —0.1, —1, respectively. 

{ Our space curve consists of three branches, making with each other angles 
of x/3 at their point of intersection, z=0, as is known from the theory of con- 
formal mapping. 

t Kempner, Téhoku Mathematical Journal, loc. cit. (On the shape of poly- 
nomial curves.) 
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These hints are sufficient to characterize qualitatively the 
space curves.for a 2"+ --- for the case when the /otal 
number of_ maxima, minima, and accessory points of inflection 
(each accessory point given a weight 2) is the largest possible.* 
Algebraically, this total number is obtained by a careful exam- 
ination of f’(z) and f’’(z) for real z. 


Fic. 18 


(3) If our total number is less than n—1 (taking multiplicities 
into account), as in z¢-+2z2?+1, (one minimum, no points of in- 
flection, and hence in particular no accessory points of inflection, 
total number =1<3), some of these free branches will in turn 
have coalesced. As an illustration of this state of affairs, we 
consider z¢+2z2?+1. Since this function is obtained from 2 
—2z?+1 by a substitution z=iz, our whole space curve is ob- 
tained by rotating Fig. 17 through 7/2 around the wu axis, thus 


* For a polynomial of degree u, remembering that each accessory point 
corresponds to one maximum and one minimum, we have M,+M;+2A Sn-1, 
where M,=number of real maxima, M;=number of real minima, A = number 
of accessory points. The sign of inequality holds when the curve does not pos- 
sess the maximum number n—1 of extremes and some of the accessory points 
of inflection have become lost. 
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changing the hyperbola in the xy plane to one intersecting the 
y axis. 

If a continuous transformation of 24—2z?+1 into 24+ 22?+1 
is desired which shall illustrate in typical fashion how handles 
attached to the real curve z*—2z?+-1 first detach themselves 
from the real curve and are afterwards recaptured by it, a 
suitably chosen variation of coefficients passing through z‘—z? 
is preferable to our chain z*+2?, 2=>¢e= —2. In the general case, 
the situation is governed by the curve v(x, y) =0; any point in 
the xy plane where two branches of v(x, y) =0 intersect (double 
points of v=0) is the projection of a point on our space 
curve where two handles or detached branches are joined, and 
conversely. This agrees with the theory of conformal mapping 
and its breakdown at points for which f’(z) =0. At all such points 
on the space curve where handles or branches are joined (in- 
cluding the real curve), they are of the character of intersections 
of a surface at a saddle point. If a point starts moving along 
any branch of the space curve at u= +, it is always possible 
to travel along the space curve to u= — ~. (This is a conse- 
quence of the last remark combined with the Fundamental 
Theorem.) 

I hope that enough has been said about these space curves to 
justify the importance I attach to them for the problem of com- 
plex roots. A collection of strong wire models of such curves, and 
their introduction into advanced elementary work, is very de- 
sirable. 


THE UNIVERSITY OF COLORADO 
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A THEOREM ON HIGHER CONGRUENCES* 
BY LEONARD CARLITZ 


1. Introduction. Let D=D(x, p”) denote the totality of poly- 
nomials in an indeterminate x with coefficients in a Galois field 
GF(p") of order p*. Consider the congruence 
(1) i" —t=A (mod P), 
where A and P are in 9, and P is irreducible of degree k, say. 
The sum 

is congruent (mod P) to a quantity in GF(p"); we denote this 
residue by p(A). It is easily seen that the congruence (1) is solva- 


ble in D if and only if p(A)=0. A better condition is furnished 
by the following theorem. 
THEOREM. If we put 
(2) P= q2"'4+---+ 
P’ = + — +--+ + 


where c; is in GF(p"), then the congruence (1) is solvable in D 
if and only if P’A is congruent (mod P) to a polynomial of degree 
<k—1. More generally, if 


= box*® 34+ --- + (mod P), (b; in GF(p")), 
then p(A) 
In this note we give a new and direct proof of this theorem.t 


2. Proof of the Theorem. For arbitrary A (mod P) we con- 
struct the polynomial 


= — A)¢ — (mod P), 


* Presented to the Society, April 19, 1935, under a different title. 
t See L. Carlitz, On certain functions connected with polynomials in a Galois 
field, Duke Mathematical Journal, vol. 1 (1935), p. 164. 


1935-] THEOREM ON HIGHER CONGRUENCES 845 


in which the coefficient of t*—! is evidently —p(A). For our pur- 
poses it will be convenient to make use of an alternative defini- 
tion of f(t). Let x denote a root of P=0; then x defines the 
GF(p"*). Then A =A(x) is an element of the enlarged Galois 
field; f(t) is evidently the unique polynomial, with leading coeffi- 
cient =1, having the roots A™. Clearly all the coefficients of 
f(t) lie in GF(p*). To calculate them we proceed as follows. Let 


(3) = (mod P), 
i=0 
where a;;, (t,7=0,---,k—1), are in GF(p*). But the equations 


(3) evidently imply the following representation of f(t) as a de- 
terminant: 


f(t) = (— 1)*| — 


so that by the remark at the beginning of this section 
k-1 
(4) p(A) = 
i=0 
On the other hand, making use of (2) and (3), 


k 
= (co = 1), 


k k-1 
i=1 j=0 
so that the coefficient of x*~! is 
k 
(5) bo = 


Note next that (3) implies 


k-1 
= 2, (mod P) 
i=0 


k—2 k-1 


i=0 j=0 


from which it follows that 


— 
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(6) = Gi,j-1 — 
Put <=j—1, and (6) becomes 

— = G=1,---,&—1). 
Substituting into the right member of (5), we see that 


k-1 
bo = — + 
j=1 


Goo + G11 + °° + 


If we compare with equation (4), we have at once p(A) =do. 
This completes the proof of the generalized form of the theorem. 
In particular, if P’A is congruent (mod P) to a polynomial of 
degree <k—1, then b)=0, and the congruence (1) is solvable. 

3. Concluding Remark. The coefficients of f(t) are, but for 
sign, the elementary symmetric functions of the quantities 4” 
(mod P). As we have seen above, the coefficient of ¢*~! is inti- 
mately connected with the congruence (1). Similarly, the last 
coefficient 


A 1+ = { (mod P) 


is connected with the congruence 
"-l=A (mod P). 


Indeed, the method of §2 leads very naturally to F. K. Schmidt's 
proof of the theorem of reciprocity: 


where P and Q are primary irreducible of degree k and /, re- 
spectively.* 

The question arises whether the remaining coefficients in f(t) 
are connected in any direct manner with criteria for the solva- 
bility of higher congruences. 


DvuKeE UNIVERSITY 


* F. K. Schmidt, Sitzungsberichte der Physikalischmedizinischen Societat 
zu Erlangen, vol. 58-59 (1928), pp. 159-172. 
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THE HILBERT SPACE POSTULATES— 
A FURTHER REDUCTION 


BY A. E. TAYLOR 


In a recent issue of this Bulletin the writer gave a set of 
postulates for abstract Hilbert space in which equality was 
treated postulationally. The result there obtained was a set of 
eleven postulates, independent among themselves, defining a 
linear vector space of a special sort. Three more postulates were 
then added, to define a Hilbert space.* As a result of investiga- 
tions carried on by I. E. Highberg and myself concerning postu- 
late systems for general vector spaces,f it is now possible to 
make a further simplification of the Hilbert space postulates. 
The equality relation as an undefined idea of the system is not 
independent, and may be defined in terms of the other notions 
in the system. When this is done the number of postulates is re- 
duced by three. The postulate system for Hilbert space is then 
the following, where A denotes either the complex number sys- 
tem C, or the real number system R. 

1. The class K contains at least one element. 


2. Ifx, ye K, thenx+yeK. 

3. IfaeA and xe K, thena-xeK. 

4. If x, ye K, then (x, y) €A. 

5. If x, y, z€K, then (x+y, z) =(x, z)+(y, 2). 

6. If x, yeK, then (x, y) =(y, x). 

7. If aeA and x, yeK, then (a-x, y) =a(x, y). 

8. If xeK, then (x, x) =0. 

9. For each integer n>O there exist elements x1,---,*x,€K 
af and only if a,=a2= --- =a,=0. 


10. K is separable according to the norm (x, x)¥? =|]. 
11. K is complete according to the norm (x, x)!?=||x!|. 


DEFINITION. If x, yeK, we say that x is equal to y and write 
x=y if and only if (x+—1-y, x+—1-y)=0. 


* A. E. Taylor, this Bulletin, vol. 41 (1935), pp. 439-448. 
+ A. E. Taylor and I. E. Highberg, On postulate systems for normed vector 
spaces. 
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In order to prove that these postulates do define a Hilbert 
space we must establish the following propositions, which were 
postulates in the original set. 

I. If x, ye K, and tf for each ue K, (x+—1-y)+u=u, then 
x=y. 
II. If x eK and (x, x) =0, then x+y=y for each ye K. 
Ill. If x, y, we K and x=y, then (x, u) =(y, u). 
ProoF oF I. Let 


(e+ — 1-9) + + — 1-0, (@ + — 1-9) + + — = 


Then, since u is arbitrary, we may replace it by 0-4. On ex- 
panding, with the aid of postulates 5, 6, 7, we easily find that 
(x+—1-y, x+—1-y) =0, or x=y. 

Proor oF II. Let (x, x)=0. Then, just as in the previous 
proof, we obtain the equation 


+ y) + — (2 + y) + — 1-y) = 0. 


This means, by definition, x +y=y. 

Proor oF III. Let (x+--1-y, x+—1-y) =0, and let u be an 
arbitrary element of K. Then (x, u) —(y, u) =(x+—1-y, u). But 
by the Schwarz inequality, which does not depend for its proof 
on the use of proposition III, we have 


O<|(x+—1-y,u)| —1-y,a+ = 0. 
Therefore (x, u)=(y, u). 

Postulates 1-8 are independent among themselves, as is 
proved by the examples given in the previous paper. The inde- 
pendence of the set taken as a whole is an open question. 


CALIFORNIA INSTITUTE OF TECHNOLOGY 
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A NOTE ON LIPSCHITZ CLASSES 
BY E. S. QUADE 


This note consists in the application of some results of Hardy 
and Littlewood* on fractional integrals to a theorem of Paley 
and Zygmund{ and gives a generalization of that theorem. 

We consider only functions of the Fourier power series type. 
That is, f(x) is periodic in 27, integrable, and with a Fourier 
series of the form 


n=0 
In dealing with functions of the class Lip(a) or Lip(a, p), a1, 
this restriction is a matter of convenience rather than one of 
necessity. { 
A function f(x) is said to belong to the class Lip(a), where 
0<a 3X1, in the interval (—7, 7), if 


f(x + h) — f(x — h) = O(h*) 
uniformly for —rSx—h<x+hz, and to Lip(a, p), where 
p21, 0SaKX1, in (—7, 7), if f(x)eL, and 
f | f(x + h) — f(x — h) = 


Tl.e functions ¢,(é), (n=0, 1, 2,---), are the Rademacher 
functions. § 


* Hardy and Littlewood, Some properties of fractional integrals I, Mathe- 
matische Zeitschrift, vol. 27 (1927-28), pp. 565-606. We will refer to this 
paper as (HL), 

t Paley and Zygmund, On some series of functions, Proceedings Cambridge 
Philosophical Society, vol. 26 (1930), pp. 337-357. A. Zygmund, Trigonometri- 
cal Series, 1935, §5.61. We will refer to this book as (Z). It contains extensive 
bibliographical references. 

t Hardy and Littlewood, A convergence criterion for Fourier series, Mathe- 
matische Zeitschrift, vol. 28 (1928), pp. 612-634, in particular, §2 and §3.5. 
See also (Z), §7.4. 

§ For definitions and properties see (Z), §§1.32 and 5.5 to 5.61. 
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THEOREM 1. Let Co, C2, +, Cn, be a sequence of real 
or complex numbers such that 


> (log n)'+e 
n=2 


converges for an €>0. Then, for almost all values of t, the function 
n=0 
belongs to the class Lip(a), (Q0Sa<1). The theorem is false for 


the case a=0 or a=1 tf €=0. 


As a consequence of the theorem of Paley and Zygmund men- 
tioned above it follows that 


n=1 


for almost all values of ¢, is a continuous function (since the 
series converges uniformly we put f; (x) equal to the sum of the 
series). That is, we have 


i (x)e Lip (0). 
If by the symbol f;-.(x) we denote the integral of ff (x) of order 
a, we have* 

fi.a(x)e Lip (a). 
But 


=i >, 
f(x) — Copo(t). 


To show that the theorem is not true in the case a=1, for 
€=0, we consider the function 


* (Z), §§9.80 and 9.81. A function satisfies a condition Lip*(«) or Lip*(a, p) 
when it satisfies a condition analogous to that for Lip(a) or Lip(«, p) but with 
o small in place of O large. In each of our theorems Lip(a) or Lip(a, p) may 
be replaced by Lip*(a) or Lip*(a, p), respectively, except in the case a=1; 
this follows from Theorems 18, 21, and 22 of (HL). 
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x 2™m log (m + 1) 


This can not belong to Lip(1) for any sequence of signs since it 
is the integral of 


+ 


m log (m + 1)’ 


which is Paley and Zygmund’s example* of a series which does 
not represent a bounded function for any sequence of signs. 

For the case of Lip (a, p), (OSaS 1/p= 1), we have a similar 
theorem. 


THEOREM 2. Let Co, C2, , Cn, be sequence of real or 
complex numbers such that cy converges. Then, for al- 
most all values of t, the function 


f(x) 


n=0 


belongs to the class Lip (a, p), (p21, OSaS1). 
Since | is convergent, it follows thatT 


fi (x) 


n=1 


belongs to Ly, (p21). 
Now by a theorem of Hardy and Littlewoodt 


Lip (a, p). 
But, as before, 
fea(*) = fulx) — cobo(t). 


As a corollary of the following theorem we have a better 
theorem for the case 1<pS2. 


* Paley and Zygmund, loc. cit., p. 350. This gives the case a=0, e=0. 
+ (Z), §5.6 (iii). 
t (HL), Theorems 21, 22, and ff. 
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THEOREM 3. If 
+2 
> (1 < p’ 2), 
converges, then 


+00 
f(x) encima 


belongs to the class Lip (a, p), 1/p+1/p’=1, (OSaZ1). 


From the Young-Hausdorff* theorem we have 
(Lf 
Td n=—2» 


+00 
f(x) 


n=—2 


where 


Since f*(x)eL,, we havef 
+c 
— = f#(x) = € Lip (a, p), (p2 2). 


Coro.iary. If 


converges, then 


+20 
f(x) >. 


belongs to Lip (a, p) for every p such that 1S pS2. 


This follows because 


1 1/p 1 1/2 


Brown UNIVERSITY 


* (Z), §9.1. 
+t (HL), Theorem 21. 


n=—20 
n=—0o 
+00 
> , 1), 
n=—20 
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A NEW PROOF OF THE EQUIVALENCE 
OF E. H. MOORE* 


BY J. P. BALLANTINE 


The theorem which is here called the equivalence of E. H. 
Moore is a statement of a necessary and sufficient condition 
for the existence of a solution of a system of linear equations. 

In Moore’s own proof of this equivalence, there is a hypothe- 
sis to the effect that if the Hermitian square of a vector a’ is 
zero, all the elements are zero. That is, 


ava, = = 0’. 


In plain language, take a set of complex numbers, or numbers 
in any number system, multiply each by its conjugate, and add; 
then if the sum is zero, all the original numbers are zero. 

This hypothesis is satisfied in many number systems, and 
holds in every case in which Moore used his equivalence. Un- 
fortunately, the writerf used the equivalence in the number 
system of residues modulo 2. For such a number system the 
hypothesis does not hold, as for example in the case of the vector 
a’=(1, 1), where a’a,=1+1=0 (modulo 2). To form a logical 
basis for this work, it was necessary to have a proof, such as the 
following, of the equivalence which would not use the limiting 
hypothesis. 

A general system of linear equations is written 


ie 


in matrix form. The prime range (’) and the second range (’’) 
are conceptually distinct, so that the matrix a/’ may be square 
or not. Let ” be the number of rows in the matrix, or the number 
of equations in the system, and m the number of columns, or 
unknowns. 

The symbol 3 is read “there exists,” and 35 is read “such 
that.” The statement that a solution of the system of equations 
exists is written 


* Presented to the Society, August 30, 1932. 
+ A postulational introduction to the four color problem, Six Studies in Mathe- 
matics, University of Washington Publications, 1930, vol. 2, No. 1, p. 9. 
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= (X) 
Two vectors a’ and 8, are said to be orthogonal if 
a’B, = 0. 
A vector a’ is orthogonal to all the columns of a matrix B/’ if 
= 0". 
It may be shown by straightforward substitution that 
where Y is the proposition 
= 0"->-7'n, = 0. 


That X implies Y and, conversely, Y implies X, constitutes 
the equivalence of E. H. Moore.* In Moore’s proof of the con- 
verse, he normalizes the columns of the matrix, and for this step 
he needs the hypothesis referred to earlier. The following proof 
does not require this hypothesis. 

To prove that Y implies X, assume Y true and X false. Then 
the proof is complete if we are led to a contradiction. 

Commutative multiplication is assumed here, but the proof 
can easily be extended to a noncommutative number system. 
Unique division, except by zero, is also assumed. 

A set, is built up of vectors Y1,m, Y2,m; °° » Which, 
together with all their linear combinations, will include all possi- 
ble vectors orthogonal to all the columns of a/’. There appears 
to be no obvious way to write down the set I, at one step, but 
it can be done by mathematical induction. Two sequences of 
sets of vectors are formed 


A; is the set of vectors consisting of the first k columns of a/’, 
and I; is a set of vectors which, together with all their linear 


combinations, includes all the vectors orthogonal to all the vec- 
tors of Ax. 


*E. H. Moore, Lectures on General Analysis. Memoirs of the American 
Philosophical Society, vol. 1, 1935. 
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Since Ao is a null set, the vectors of ') may be written down, 
by inspection, as the rows of Hilbert’s unit matrix 6’, namely, 


where 


= 1,0,0,---,0, = 0,1,0,---,0, 


The formation of I'¢.41) from T; is as follows. All the vectors 
of I’, are assumed orthogonal to all those of Ax. Now A «41 in- 
volves one new vector v, not in Ax. To see if all the vectors of 
I’, are accidentally orthogonal also to v,, we consider the follow- 
ing NV; numbers: 


0; = v5.0", (j= 1,2,---,M,), 


where WN; is the number of vectors 71,4, --- , Yv.xin 

Case 1. Now if all the numbers Q; are zero, then all the vec- 
tors of I, are orthogonal to all the vectors of A (41). This will 
happen, of course, if v, the (k+1)st column of a/’ is linearly 
dependent on the preceding columns. It may, as far as we have 
yet shown, happen in other cases. 

CasE 2. All the numbers Q; are not zero, so that Q; +0, where 
J is some particular value of 7. Let P denote the reciprocal of 
Q,;. In our number system such reciprocals are supposed to exist 
uniquely. As a tentative definition of T'.41, consider 


To test that all the vectors },4.41) are orthogonal to all the 
vectors of A(z+1, note that it is obvious for the first k of the 
latter. In the case of v,, 


Q; — = 0, 


by the definition of P. In testing the linear independence of the 
vectors of 441) as defined above, we find that one of them, 
3.(h41), is identically zero. That one must be deleted, and the 
others constitute the set I'(z,:) as finally defined, so that the 
number of vectors in T'¢,.41) is one less than in ’;, or 


= 
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N (e+1) = Ni 


The linear independence of the N;,—1 vectors of T'x+1) fol- 
lows directly from the linear independence of the vectors of T;. 
Thus the rank of I; is N;, the same as the number of vectors 
in it. Let the rank of A; be R;. Then in Case 1, 


Nut» = Ni, = Ri, or Re +1, 


according as rv, is or is not linearly dependent on the vectors of 
A,. In Case 2, 


Nosy = 1, = Re +1. 
Therefore, in both cases, 
N + 2 Ne t+ Re. 


When k=0, No=u=the number of rows in a/’ and Ro=0; 
so that No+Ro=2, and therefore N.+R,=n. 

Now prolong the sequence Ao, - - - , Am, by defining A (m41) as 
the set of all columns of a/’ and the vector y,. Then I'm41) is 
formed from [,, in the manner described for forming Tas.) 
from 

To obtain the contradiction, we are assuming Y true and X 
false. Since Y is true, all the vectors of I, are orthogonal to all 
the vectors of Aim, so we have Case 1, where I'ni1)=I'm. 
Therefore N m+ = pt 

Since X is false, the last vector of A(m41) is not expressible 
linearly in terms of the vectors of Am. Therefore Romi1) =Rat+1, 
so that 


N + > Na + Rn 


Now the sets A (m41) and (m4) are mutually orthogonal, and 
the vectors of those sets have only 2 elements. It is, therefore, 
impossible for the sum of their ranks to exceed n, and we have 
the desired contradiction. 

Or, by extending the sequence of A’s further, one could ob- 
tain as a contradiction a single set of vectors in m-space with 
rank +1. This, of course, could be reduced step by step to a 
set of vectors in one dimension of rank two. 


THE UNIVERSITY OF WASHINGTON 
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ON THE BERNOULLI DISTRIBUTION* 


BY SOLOMON KULLBACK 


1. Introduction. The Bernoulli or binomial distribution, of 
central importance in the theory of mathematical probability 
and statistics, has been the subject of considerable study. The 
derivation of the moments, or recursion formulas for the mo- 
ments, of this distribution, as usually presented, involve the 
use of moment-generating or characteristic functions, or the ex- 
plicit form of the distribution itself.t In the following we will 
derive these moments in an elementary manner and extend the 
results to the Poisson exponential distribution, distributions of 
the Lexis and Poisson types, and the multinomial distribution. 


2. Preliminary Notions. We need and use the following as- 
sumption, definitions, and theorem. 

EMPIRICAL AssuMPTION. If an event which can happen in two 
different ways be repeated a great number of times under the 
same essential conditions, the ratio of the number of times that 
it happens in one way, to the total number of trials, will ap- 
proach a definite limit, as the latter number increases indefi- 
nitely.f 

DEFINITION. The limit described in the empirical assumption 
shall be called the probability that the event shall happen in 
the first way under those conditions.§ We shall express the fact 


* Presented to the Society, November 30, 1934. 

t See, for example, A. Fisher, The Mathematical Theory of Probabilities, 
2d ed., p. 104 ff.; H. L. Rietz, Mathematical Statistics, 1927, p. 26 ff.; V. Mises, 
Wahrscheinlichkeitsrechnung, 1931, pp. 131-133; Risser and Traynard, Les 
Principes de la Statistique Mathématique, 1933, pp. 39-40 and 320-321; V. 
Romanovsky, Note on the moments of the binomial (q+p)" about its mean, 
Biometrika, vol. 15 (1923), pp. 410-412; A. T. Craig, Note on the moments of a 
Bernoulli distribution, this Bulletin, vol. 40 (1934), pp. 262-264; A. R. Cra- 
thorne, Moments de la binomiale par rapport al’ origine, Comptes Rendus, vol. 
198 (1934), p. 1202; A. A. Krisknasuami Aygangar, Note on the recurrence 
formulae for the moments of the point binomial, Biometrika, vol. 26 (1934), pp. 
262-264. 

tJ. L. Coolidge, An Introduction to Mathematical Probability, 1925, p. 4. 

§ J. L. Coolidge, op. cit., p. 4. 
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that the relative frequency f/s (the ratio of the number of times 
the event happens in one way to the total number of trials) 
is assumed to have the probability » as a limit when s—, in 
abbreviated form by writing E(f/s)=p, where E(f/s) is read 
“expected value of f/s.”* 

THEOREM OF COMPOUND PROBABILITY. If a compound event 
consists in the conjunction of any number of independent 
events, the probability of the compound event is the product 
of the probabilities for the individual events. f 


3. Bernoulli Distribution. Suppose that in m independent 
trials we observe an event to occur x times. Then, applying 
the notions of §2, we have 


x x(x—1) 


a(x—1)--- 
n(n—1)--- (n—r+1) 


by 


by 


n—x)(n—x—1) -- - (n—x—r+1) 

n(n—1)--- (n—r+1) ) 
x(x—1) - - - (w—r+1)(n—x)(n—x—1) - - - (n—x—s+1) 

n(n—1)(n—2) -- - (n—r—s+1) ) 


E 


= rr 
for r+s<n. Each fraction in (1) represents the ratio of observed 


r-plets to the total possible number of such 7-plets. We may also 
write (1) as 


(2) {E(x/n)¥ {E((n — x)/n)}" = E([x] [n — ), 


where [x]’ is the factorial function} x(«—1)(x—2) - - - (x—r+1) 


* H. L. Rietz, op. cit., p. 9. 

7 J. L. Coolidge, op. cit., p. 18. 

t Whittaker and Robinson, The Calculus of Observations, 1924, p. 8. See also 
J. F. Steffensen, Factorial moments and discontinuous frequency functions, 
Skandinavisk Aktuarietidsskrift, vol. 6 (1923), pp. 73-89. 
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and {E(x/n)}" is the rth power of the expected value of x/n, 
and r+s <n. 

In order to find the expected value of the rth power of the 
variable (or the rth moment about the origin), we avail our- 
selves of the fact that 


(3) a,.o[x] + + des, 
where a;,;, (i, 7=0, 1, 2,--- ), is the element in the ith row and 
jth column of the table* 
0 0 0 0 0 0--- 
1 1 1 
3 7 15 31 63--- 


nN 
on 
‘Oo 


in which @0,,=0, ((=0, 1,---), aio=1, 
(t=1,2,---). We thus have 


E(x") = ay, 


+--+ + 41,rE(x), 


(4) 


or 
(5) E(x’) = a;,o[n] pr + + 


For the expected value of the rth power of the deviations of x 
from its expected value (or the rth moment about the expected 
value), we have 


(6) E{(« — E(x))"} = (— 1)°C,, E(x) 


in which we replace E(x’) by its value as given in (5). This 
yields 


* Whittaker and Robinson, op. cit., p. 9. 


_ = 
} 

} 
| 
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(7) Ef{(x — E(x))"} = + + Bip, 


where 


B, = (-— 
i=0 
r—j—-1 
B,-; = >> (— 
i=0 
4. Poisson Exponential Distribution. The moments of the 
Poisson exponential distribution are derived from those of 
the Bernoulli distribution by letting n—« with lim,..E(x) 
=lim,..1p=m, where m is finite. Accordingly, we have from 
(4) and (7), respectively, 


(8) E(x") = a;,om? + + + 
(9) Ef{(x — E(x))"} = + A,ym™ +--+ + Am, 


where 


r r—F—1 
i=0 t=0 


5. Lexis Type Distribution.* Consider N independent sets of 
n independent trials each, with the probability for the occur- 
rence of an event constant for a set, but varying from set to set. 
Let the varying probability be fi, po, - - - , py, and the number 
of occurrences observed, respectively, x1, %2,-°-- , xy. 

Then, in a manner similar to that of §3, we have 


(10) E([x]” /[nJ) = y where [x]” = 
i=1 


If we set p =)" _b/N, we have, corresponding to (5), 
E(x?) = p + + --- 
+ 
where x = x?/N, and corresponding to (7), 
(12) E{ (x® E(«%™))r} = Bp) + +--+ + Byp™, 


Sincet 


*H. L. Rietz, op. cit., Chap. 6. 
7 Chrystal, Textbook of Algebra, 2nd ed., vol. 2, p. 49. 


1 N rT 
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we see, by comparing (5) and (11), that 


(13) E(x") 2 Ex(x"), 


where E; represents the expected value from a Lexis type dis- 
tribution and Eg the corresponding expected value from a 
Bernoulli distribution with _pi/N. Indeed E1{(x+a)"} 
> Ex{(x+a)} for a=0. 


6. Poisson Type Distribution.* Consider a set of n inde- 
pendent trials with the probability for the occurrence of 
an event varying from trial to trial. Let these probabilities be 
Pi, po, + - - , Pn. The study of the moments for this case is really 
the consideration of the moments for the distribution of the 
total of a Lexis type distribution of m independent sets of one 
trial each. Thus, if the number of observed occurrences of an 
event in x trials is x=x,+x2.+ --- +x,, where x;=zero or one, 
then 


(14) E(x) = pi E(x;) = > Pi, 


(15) E(x") = pit + a 


where a,” ++ $x!) for sitse+ ---+s,=r, 
(t=1, 2,---, Rk). The right member of (15) follows from the 
fact that E(x,;7)=;, as may be readily seen from (4) for n=1 
and r22. 

For the moments about tke expected value, we have 


r 


(16) E{(x — E(x))"} = (— 19°C, { E(x) }*, 
t=0 
where E(x’) is given by (15) and { E(x) }"=(pitp.+ 
If we repeat the preceding argument with constant probabil- 
ity p=).”_,bi/n, we find that 


(17) E(x) = a, np + a2C, 2p? + + 


* H. L. Rietz, op. cit., Chap. 6. 
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Since* 
1 n r 
N 
we see by comparing (17) and (15) that 
(18) Ep(x") = Ep(x’), 


where Ep represents the expected value from a Poisson type 
distribution and Eg the corresponding expected value from a 
Bernoulli distribution with Indeed, (x+a)"} 
>Ep{(x+a)"} for a=0. 

A comparison of (17) and (5) shows that 


(r) k (r) 
(19) C105 = Qk,r—k» or = Riaz 


7. The Multinomial Distribution. Consider a trial in which 
one of r mutually exclusive events may occur. Let the re- 


spective probabilities of occurrence be f1, p2,---, p, with 
Pitpet+ --- +p,=1. If in nm independent trials there are ob- 
served x1, X2,- occurrences, respectively, of the 7 events, 
then 


(20) E{ = pipe 


Thus, for example, E(x,;) =n(n—1)p:p; and E(x;x;) — E(x;)E(x;) 
=n(n—1)p:p; -—npnp;= 
Now we may writet 


xt => 


A=0 


where A>0! are the differences of zero and 


1 
AO! = — — 1)? + —1)(A—2)'—--- 
* This can be shown by using the fact that f(x1, x2, , Xn) =) 
subject to the condition --- +x, =constant, has its maximum value 
for x) =xX2= +--+ 


7 J. F. Steffensen, Interpolation, 1927, p. 13, equation (3). 


= 
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Accordingly, 
(2 1 ) 


[ x, > 


«=0 K ! A=0 ! ! 


and 


(22) E(x1x2 ) 


K+A+-++t+u 


AO... A’0" [n] Pr 


ll 


where A*0;'=1 for 14+7 and A*0;'=A*0". The latter result is 
true for any m and k, l,---, m since [n|*+---+#=0 for 
K+A+ --- and A*0"=0 for a>r. 

The procedure with respect to the product moments about 
any constant is similar to that outlined. Thus, since* (x;—a;)! 


=) [xi we have 
(23) (%1 — a) "(a2 — a2)! -- (x, — a,)™ 


-Ae(—a,) [xi] [xe] - 


and 
E{ (x1 — a1)*(x2 — a2)! (a, — ay)™} 
(24) 
= (1 + + pods + --- + p,A,)” 
a:)"(— -- (— a,)*, 
where di, are any constants and A#(—a;)'=1 for 


i~j and A# (—a;)'=A*(—a;)". Sincef A=e? —1, where D=d/dx, 


* J. F. Steffensen, op. cit., p. 13. 
¢ Whittaker and Robinson, op. cit., p. 63. 


— 
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(25) E{ (x — — @s)'--- (% — a,)™} 


= (pre t+ pe + pe 


where D,=0/0x;. 


Wasuincrton, D. C. 


ON A RESULTANT CONNECTED WITH 
FERMAT’S LAST THEOREM 


BY EMMA LEHMER 


E. Wendt* seems to have been the first to introduce the re- 
sultant of x*=1 and (x+1)"*=1 in connection with Fermat’s 
Last Theorem. This resultant can be expressed by means of the 
following circulant of binomial coefficients 


1 aks Cons 
A, = 


In his book on Fermat’s Last Theorem Bachmannj{ proved 
that if p is an odd prime and if A,_; is not divisible by p*, then 
Fermat’s equation x?+y”?+2? =0 has no solution (x, y, 2) prime 
to p. 

S. Lubelskyf proved in a recent paper, using the distribution 
of quadratic residues, that if p27, A,-1 is not only divisible 
by #*, but by p%, thus annulling Bachmann’s criterion except for 
p=3 and p=5. 

We shall now show how, by a straightforward manipulation 
with the above determinant, one can prove much more. 


THEOREM 1. A,_; is divisible by p?~*q2 for every prime p, where 
ts the Fermat quotient (2?-!—1)/p. 


* Journal fiir Mathematik, vol. 113 (1894), pp. 335-347. 
¢ Das Fermatproblem, 1919, p. 59. 
} Prace Matematyczno-Fizyczne, vol. 42 (1935), pp. 11-44. 
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Proor. First add to each element of the last column of A,_; 
the corresponding element of all the other columns. The ele- 
ments of the last column now become equal to 


1+ + Cota = — 1. 


Next increase each element of the first —3 columns by the 
element immediately to the right of it. All the elements of the 
first p—3 columns are now of the form 


+ = = pli, (k =0,1,---,p— 2). 


Since, as is well known, J; is an integer for p a prime, it follows 
that » is a factor of each of the first —3 columns. Also 
(2?-1—1) comes out of the last column, and hence A,_, is di- 
visible by (2?-!—1)p?-*=p?-*q2, which is the theorem. 

For example, if p=5, 


5 10 6 15 
5 4 15 
10 10 4 15 


It is interesting to notice that although Theorem 1 is in no 
way dependent on the solvability of Fermat’s equation, never- 
theless it enables us to replace Bachmann’s criterion by the fol- 
lowing one. 

If A,-1 is not divisible by p?—, then x?+y?+2?=0 has no solu- 
tion (x, y, 2) prime to p. 

This in fact is merely a restatement of Wiefrich’s criterion,* 
which states that if p does not divide g2, Fermat’s equation has 
no solution in integers prime to p. 

For example, since we have seen that A, is not divisible by 54, 
it follows from our criterion that x5+-y5+25=0 has no solutions 
(x, y, 2) prime to 5. 


* Journal fiir Mathematik, vol. 136 (1909), pp. 293-302. 
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Before proceeding further it may be of interest to give a short, 
though perhaps less elementary proof of Theorem 1, based di- 
rectly on the definition of the resultant of two polynomials as 
the product f(a:)-f(a@z) - - - f(an), where f is one polynomial and 
the a’s are roots of the other. 


p—l 
Apa = + 1], 
j=1 
where ¢; are all the (—1)st roots of unity. 
For €;=1, we get the factor 2?-!—1, and for e,= —1, we get 
—1. For ¢;=€, any complex root of unity, we have 


(e + C5--1,1€ Cp-1,2€7 
But since, for p a prime, 
= (— 1)* + 


where the c’s are integers,* we have 


(e+ 
= ee? — 1)/(e+ 1) + = pf), 


where f(x) is a polynomial with integral coefficients. Hence 


Apa = — — 1) = — 1)-p?]] fled), 


where |] f(e;) is an integral symmetric function of the roots of 
(x?-!—1)/(x?—1) and hence an integer. Hence A,_; is divisible 
by p?~7q2. 

In comparing Theorem 1 with that of Lubelsky we see that 
Theorem 1 says more, except for =7, in which case Theorem 1 
guarantees divisibility of Ag by 7° instead of 78. On examining 
this case more closely we find that as a matter of fact Ag=0. 
Indeed, we have in general the following theorem. 


THEOREM 2. A, =0 if and only if n=6k. 


Proor. In order that A, =0 it is both necessary and sufficient 
that x"=1 and (x+1)"=1 have a root p in common. But the 
roots of x"=1 are the mth roots of unity. Hence we can write 


* Lucas, American Journal of Mathematics, vol. 1 (1878), pp. 229-230. 
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p=cos 6+7 sin 0, but since at the same time p is a root of 
(x+1)"=1, we have 


|p +1| = 1 = (cos@+ 1)? + sin? = 2+ 2 cos@, 


or cos = —1/2 and @= +27/3. This condition will be satisfied 
if and only if p=w or w’, while (9 +1) = —w? or —w, and hence 
(o+1)"=1 if and only if m is a multiple of 6. 

One can easily show by adding to each element of the first 
and second column of A¢, (written in determinant form) the 
corresponding elements of every third column, that the elements 
of the two” resulting columns} will, be equal and hence that 
Ac. =0, but the writer has not been’able to show from the circu- 
lant definition of A, that A, 0 if 1 is not a multiple of 6. 

In conclusion we give a short table of A,_1. 


p Ay-1 


—3=-—(2?-1) 

5  —375= —(24—1)-5? 

7 0 

11 —210 736 858 987 743 = —(2!°—1)-118-31? 

13 0 

17 —1 562 716 604 038 367 719 196 682 456 673 375= 

— (216—1)-1714. (33-5- 73-257)? 
19 0 
It appears from this table of A,_1, and can be shown without 

difficulty for any even n, that A, is —(2"—1) times a perfect 
square. It can also be shown that Aq divides A, if d divides n. 


BETHLEHEM, PENNSYLVANIA 
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CORRESPONDENCES CONNECTED WITH 
A PENCIL OF »n-ICS* 


BY A. R. WILLIAMS 


1. Introduction. All curves of order m passing through 
(n/2)(n+3)—1 points pass through (m—1)(m—2)/2 other 
points. That is, the m? base points of a pencil of m-ics are de- 
termined’ by the former number. When (/2)(n+3)—2 are 
fixed in general position, and another moves on a curve of 
order m, the locus of the remaining (n—1)(m—2)/2 is a curve 
of order m(n?—1), which has a multiple point of order mn at 
each of the fixed points.t In a previous paper I have consid- 
ered the case when a number of base points are fixed and the 
others necessary to determine a pencil are taken consecutive 
on some curve.f{ One other situation is perhaps worth brief no- 
tice. Of the nq intersections of two curves of order m and gq, 
n>q, nqg—(q—1)(q—2)/2 can be taken arbitrarily, the others 
being then determined. Moreover, 


nin t+ 3)_, 


2 2 2 


This means that if mg—(q—1)(q—2)/2 base points are taken 
on a curve of order g, and (1/2)(n—q)(n—q+3)—1 others 
not on the latter, then if one other, say P, describes some 
locus, (¢—1)(q—2)/2 will be fixed on the curve of order gq, 
and the remaining (1/2)(m—1)(m—2)—(1/2)(q—1)(q—2) 
=(1/2)(n—q)(n+q-—3), variable with P, will, for any posi- 
tion of P, lie on the curve of order n—g determined by P and 
the second set of fixed points just mentioned. It is the purpose 
of this paper to discuss the locus of the variable points when 
P describes a given curve. 


2. Order of the Locus and Its Singularities at the Base Points. 
In this connection we may use the rational surface whose plane 


* Presented to the Society, April 6, 1935. 
+ Milinowski, Journal fiir Mathematik, vol. 77, p. 263. 
t Williams, this Bulletin, vol. 36, p. 133. 
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sections correspond to the »? m-ics that have in common the 
nq fixed points on the g-ic and all but one of the others. It will 
be convenient to refer to the set on the g-ic as the points A, 
and the others, (1/2)(m—q)(m—q+3)—2 in number, as the 
points B. The remaining fixed point, which reduces the multi- 
plicity of m-ics to a net, may be called Q. The points B and Q 
determine a pencil of (w—q)-ics. The order of the surface is 
(1/2)(n—q)(n+q—3)+2. It has a multiple point of order 
q(n —q) corresponding to the g-ic. The sections through the mul- 
tiple point correspond to the net of (n—q)-ics determined by the 
points B. In particular, to the points A on the q-ic correspond nq 
lines of the surface which pass through the multiple point. To a 
pencil of n-ics determined by the points A, B, Q, and P cor- 
responds a pencil of plane sections whose axis is Q’P’, the latter 
being the points of the surface that correspond to Q and P. The 
axis Q’P’ meets the surface in (1/2)(m—gq)(n+q—3) other 
points to which correspond in the plane the remaining base 
points of the pencil of u-ics. As P moves on a general line /, P’ 
moves on an n-ic, which is skew unless g=n—1, and has a mul- 
tiple point of order g at the multiple point of the surface. Thus 
the cone generated by Q’P’ has a q-fold edge. It meets the sur- 
face in the skew -ic and a curve whose image in the plane is the 
locus sought. To find the order of the latter we have only to 
note that to a second line /’ of the plane corresponds a similar 
n-ic on the surface which meets the cone described by Q’P’ in 
n? points, g? of which are at the multiple point of the surface, 
while another corresponds to the intersection of I] and /’. Hence 
as P moves on /, the variable points fill a curve L of order 
n? —q*?—1 which passes n —gq times through each of the points A, 
n times through each of the points B, and m times through the 
point Q, each passage through Q corresponding to an intersec- 
tion of the locus of P’ with the tangent plane to the surface at 
Q’. Thus, of the intersections of L with the q-ic, ng(m—q) are ac- 
counted for at the points A; and the remaining g[q(n—q)—1] 
correspond to the gq intersections of / with the q-ic. When 
P takes a position on the q-ic, the pencil of u-ics becomes 
the g-ic and the pencil of (n—g)-ics determined by the points 
B and Q. One of these passes through P and meets the q-ic 
in g(a—g) —1 other points which correspond to P. The others 
are the (1/2)(—q—1)(m—q—2) remaining base points of the 
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pencil of (m—g)-ics. The sum of the last two expressions is 
(1/2)("—gq)(n+q—3), the number of points variable with P. 
The remaining base points of the pencil of ( —q)-ics correspond 
to the points in which the line joining Q’ to the node of the sur- 
face meets the latter again. Each of them is a g-fold point on L, 
since P in describing / meets the g-ic q times. 


3. Self-Corresponding Points and Branch-Points. The locus of 
of points P such that one of the corresponding points coincides 
with P is the jacobian of the net of n-ics determined by the 
points A, B, and Q. It is the image of the curve of contact 
of the tangent cone drawn to the surface from Q’. This jacobian 
contains the g-ic. The residual J, of order 3(n—1)—g, meets 
l in those points which are points of L. Of the other intersec- 
tions of J and L, n—g are accounted for at each of the ng 
points A, and 2m at each of the points B and at Q. This leaves 
(n —q)(2n?+ 2ngq —6n —3q—q?—2) —2(n—3) points of intersec- 
tion which do not take place on / or at any of the base points. 
These are points at which two of the points corresponding to a 
position of P coincide. That is, they correspond to branch-points 
on /. In connection with the correspondence between the points of 
land L y.e may apply the formula of Zeuthen, 7 — n’ = 2a(p’ —1) 
—2a’'(p—1). Here a and a’ are 1 and (n—g) (n+q—3)/2, while 
n and n’ are the numbers of branch-points on / and L, and p and 
p’ are their genera. That is, 7 is the number last determined, 7’ is 
0, and pis 0. Thus p’, the genus of L, is (n/2) —a’+1, which is 
found to be (1/2) (m—q)(2n?+ 2nq—gq?—4q—7n+1) —(n—4). It 
is of some interest to verify this number. The line /, the locus of P, 
and the corresponding locus L meet in n?—q?—1—3(n—1)+ q or 
(1 —1)(m—2)—g(q—1) points which are not on J. These come 
in pairs, each pair corresponding to an apparent double point of 
the locus of P’ as seen from Q’. Moreover, with each such pair 
are associated a’ —1 or (1/2)(m—q)(n+q—3) —1 nodes on L. For 
any one of the a’+1 associated points determines the others. 
Let P; be a point of / for which one of the corresponding points, 
say R,, falls on /. When P describing / reaches R;, one corre- 
sponding point is at P; and the a’ —1 others coincide with points 
corresponding to P;, producing a’ —1 nodes. Now the genus of L 
found by considering the singularities at the points A, B, and Q, 
the (1/2)(~—q—1)(n—q-—2) q-fold points common to the pencil 
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of (n—gq)-ics, and the (1/2)[(m—1)(n—2)—q(q—1) ](a’—1) 
nodes just mentioned, agrees with the above value of p’ 
given by the formula of Zeuthen. 

If now P moves on a curve ¢ of order m, the locus of the cor- 
responding points is a curve C of order m(n?—q?—1) which 
has an mn-fold point at each of the points B and at Q, and 
an m(n—q)-fold point at each of the points A. Each of the re- 
maining base points of the pencil of (~—gq)-ics is an mg-fold 
point of C. The number of branch-points on c is also m times 
the number on /. Zeuthen’s relation then gives for the genus of 
C, p’ =(n/2) —a’+a’p+1. Here p is the genus of c, and 7, the 
number of branch-points on c, is m times the number of branch- 
points on /. This shows, as is geometrically evident, that each 
node of ¢ gives a’ nodes on C. The curves c and C meet in 
m(n* —q?—1)—m(3n—3—g) points not on J. The last number 
may be written (mn —1)(mn—2) —mgq(mq—1)—(m—1)(m—2), 
which shows that it is twice the number of apparent double 
points of the locus of P’. With each such apparent double point, 
or with each such pair of intersections of c and C, are associated 
a’—1 nodes of C. Moreover, as just remarked, each of the 
(1/2)(m—1)(m—2)—p nodes or cusps of ¢ gives a’ nodes or 
cusps on C. Then, considering the singularities of C at the points 
A, B, and Q, the mq-fold points at the remaining base points of 
the pencil of (~—gq)-ics, and the two classes of simple nodes 
just described, the genus of C is found to agree with the above. 
It is easily shown that if the locus of P passes through a point 
B or through Q, the corresponding locus C loses the n-ic of the 
net that has that point for double point. But if the locus of P 
passes through a point A on the q-ic, C loses the (w—g)-ic that 
is determined by the points B, Q, and the point A in question. 


4. Particular Cases. The cases g=n—1 and g=n—2 are of 
interest. When g =” — 2, there are only three points B. The num- 
ber of base points variable with P is 2n—5. If P describes a 
line, the order of their locus is 42—5. The order of the rational 
surface is 2n—3, and its multiple point is of order 2” —4, corre- 
sponding to the fact that the only base points of the pencil of 
conics are the three points B and Q. The conics of the net de- 
termined by the points B are the images of sections by planes 
through the multiple point. The double curve of the surface 
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consists of three (7 —2)-fold lines passing through the multiple 
point. Their images are the lines of the plane joining the points 
B in pairs. That is, through the points A and B pass a web of 
n-ics, and the m—2 remaining intersections of an n-ic of this web 
with the line joining two points B are determined by one of 
them. When P moves on the line B,Be, the locus of the variable 
points is B,B, taken n—3 times and B;Q taken m—2 times. 

When g=n-—1, there are no points B. When Q is fixed, the 

-number of points variable with P is n—2. They lie on the line 
QP; and if P moves on a general line, their locus is of order 
2(m—1), and has an n-fold point at Q and simple points at the 
points A on the (m—1)-ic. The rational surface is of order 1, 
and has a multiple point of order n—1. Sections through the 
latter correspond to the lines of the plane. If P describes a line 
passing through Q, the locus of the variable points is of course 
that line. But in this connection one fact may be noted. To 
the linc jciaing Q and A; on the (m—1)-ic corresponds on the 
surface a plane (m—1)-ic whose plane contains the line a; of the 
surface that corresponds to A;. Hence, for any position of P on 
the line QA; one of the associated points is at A,. That is, the 
pencil of m-ics have a common tangent at A;, the one having a 
double point being the (w—1)-ic and the line. 

Since the number of points A on the (m—1)-ic is just sufficient 
to determine a web of n-ics, we may take Q, P consecutive on 
some curve c of order m. That question for the base points in 
general position was considered in the article referred to at the 
beginning of the paper. The modification to be made may be 
indicated briefly. To c corresponds on the surface a curve c’, 
of order mn, having an m(n—1)-fold point at the (m—1)-fold 
point of the surface. To a pencil of plane m-ics tangent to c cor- 
responds a pencil of plane sections whose axis is tangent to c’ 
and meets the surface in n—2 other points, giving a curve on 
the surface whose image in the plane is the locus sought. The 
tangents to c’ form a developable of order 2(mn—1+), where 
pb is the genus of c and c’. For the sake of simplicity c is sup- 
posed to have no cusps and to contain none of the base points. 
Projecting c’ from the multiple point of the surface, which we 
have seen is of order n—1 for the surface and of order m(n—1) 
for c’, we find that the developable has there a multiple point 
of order 2m(n—1). The developable intersects the surface in the 
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cuspidal edge counted twice, and in a residual curve whose 
image is the locus sought. To find the order of the latter we 
have only to note that to a line of the plane corresponds on the 
surface a plane n-ic which has an (n—1)-fold point at the multi- 
ple point of the surface. Multiplying the order of the develop- 
able by n, and subtracting the 2m(n—1)? intersections at the 
multiple point and the 2m on the cuspidal edge, we find that the 
order of the locus sought is 2(2mn —2m—n-+np}). This holds for 
m=2. The locus passes 2(m+p—1) times through each of the 
points A. If the points A were in general position, the order r of 
the developable would remain the same; but the order of the 
locus would be rn—2m, that is, greater by 2m(n—1)?, and it 
would pass 2m(n—1) more times through each of the points A. 
In addition to the singularities at the base points the locus has a 
number of nodes corresponding to those intersections of the dou- 
ble curve of the developable with the rational surface that do 
not occur on the cuspidal edge. Since the genus of a plane section 
of the developable is the same as that of the cuspidal edge, 
the order of the double curve is (1/2)(r—1)(r—2) —mn— pb. If 
the cuspidal edge has a multiple point of order s, the double 
curve has there a multiple point of order 2s(s—i). There- 
fore, the double curve of the developable in question has a 
multiple point of order 2m(n—1)[m(m—1)—1] at the multi- 
ple point of the rational surface, and a multiple point of order 
4 at each of the (1/2)(m—1)(m—2) — p nodes of c’. The number 
a of stationary planes of the cuspidal edge is given by 3(r—2) 
—2mn+6p. The point of contact of such a plane is a simple 
point on both the double curve and the cuspidal edge. A point 
where the cuspidal edge is met by a tangent to itself is a cusp on 
the double curve. A curve in 3-space of order uw and genus p 
without nodes or cusps has 2(u—2)(u—3)+2p(u—6) tangents 
that meet it again. This number is reduced by 4 when an ap- 
parent double point is replaced by a node, and by the order 
uw if the apparent double point is replaced by a cusp. In the 
present instance the order of the cuspidal edge is mn, the 
genus is ~, and there is a multiple point of order m(m—1) and 
(1/2)(m—1)(m—2)—p nodes. The number y of tangents to 
the cuspidal edge that meet it again is thus found to be 
4m(1+mn—m—2n)+8+2p(mn—4). The number of inter- 
sections of the double curve of the developable with the ra- 
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tional surface that do not take place on the cuspidal edge is 
nx —2m (n—1)?[ m(n—1)—1]—a—2y—4| (1/2) (m—1) (m—2) —p], 
where x is the order of the double curve given above. Each 
such intersection corresponds to a node on the locus sought. 
Therefore, considering this number of nodes and the singulari- 
ties at the base points, the genus of the locus is found to be 
(n—3)(3mn—2m—n—3)+p(n?—8). If the base points were in 
general position, the corresponding expression would be as fol- 
lows:* (n—3)(S5mn? —6m —5n)/2+p(5n?—15n+2)/2. As is to 
be expected, both reduce to » when u =3. For we should then be 
dealing with a pencil of cubics having six fixed base points while 
O, P, the seventh and eighth, are consecutive on the curve of 
order m. Hence to each point of this curve corresponds one point 
on the locus in question. And the genus of the latter we should 
expect to be the same whether the six fixed points are on a conic 
or not. 

Finally if m==1. that is, if the two consecutive points P, Q 
describe a line 1, the loc. s of the others is that same line, and 
the above result dees not hold. But since to / there corresponds 
on the surface a rational plane n-ic having an (m—1)-fold point, 
it is seen immediately that the number of positions of P on / 
such that one of the corresponding points coincides with P is 
the number of inflections of such an m-ic, that is, 3(#—2). And 
the number of positions of P such that two corresponding points 
coincide is twice the number of double tangents, that is, 
4(n—2)(n—3). 


THE UNIVERSITY OF CALIFORNIA 


* See the paper mentioned above; 
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SOME MULTIPLICATION THEOREMS 
FOR THE NORLUND MEAN 


BY FLORENCE M. MEARS 


Absolute summability for the series }>~_,u, by the Cesaro 
mean and by the Riesz mean have been defined by Fekete* and 
by Obrechkoff,{ respectively. In each case, theorems for the 
multiplication of series summed by these means have been 
proved.{ The purpose of this paper is to establish a definition 
for absolute summability by the Nérlund mean, and to prove 
three multiplication theorems for this mean. Theorem 1 in- 
cludes Mertens’ theorem for convergent series and its extension 
for the Cesaro mean. Theorem 2 includes Cesaro’s multiplica- 
tion theorem. Theorem 3 includes the following theorem by M. 
J. Belinfante for the Cesaro mean. 

If is summable C, to U, and if vn is summable 
C, to V, and bounded C,_1, (s 20, r21), the product series > Wn 
is summable C,1, to UV.§ 

For any given series >, ,_,#x, with terms real or complex, form 
the sequence { U;}, where U; =)>oi_jttn. Let {an} be a sequence 
of positive numbers, and let A; =) ,_,@n. The series is 
said to be summable to U’ by the Nérlund mean A if 


n 
On—k+1 


lim = lim 

n> n> A; 
exists and is equal to U’.§ If , where u,, = — is 
absolutely convergent, we shall say that }>y uz is absolutely 
summable A. We shall assume that lim,....(a@,/A,) =0; then A is 
a regular method of summation.|| 


* Matematikai és Természettudom4nyi Ertesité, vol. 32 (1914), pp. 389- 
425. 

+ Comptes Rendus, vol. 185 (1928), pp. 215-217. 

t For discussion and references, see Kogbetliantz, Mémorial des Sciences 
Mathématiques, No. 51. 

§ Koninklijke Akademie te Amsterdam, Verslag, vol. 32 (1923), pp. 177- 
189. 

J Riesz, Proceedings of the London Mathematical Society, (2), vol. 22 
(1923), p. 412. 

|| Riesz, loc. cit. 
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We shall consider also the series )>,_,2%, and )>,_,w:, the 
Cauchy product of series }>*_,u, and >>,_,vx; we have the cor- 
responding sequences { V;} and { W;}. 

We shall assume that we have also a regular Nérlund mean, 
B, defined by {b,}, a sequence of positive numbers. We shall 
form the Noérlund means, defined by {cn} 
and D, defined by {d,} = 


THeoreM 1. Jf >> ,ux is summable A to U’, and in addition, 
absolutely summable A, and tf Dit is summable B to V’, then 
> ws is summable C to U'V’. 


Proor. We shall prove the theorem first with the assumption 
that V’=0. Let 


Ul =— Daw n—k+1 n—k+1) 


B, pet 


n 


n k=1 


let and be the corresponding series; 


n |n—k 
| |] vas + 
k=2 l=1 
where 
n—l+1 
p=k+1 
Hence 
< 1 J | , | 
n n—k+1 
+ ui | Ax | 1 l 
k=2 


n—1 n—k n—I+1 


} l=1 p=k+1 


n—1 
= P,|uf| + DSL| | One] + | wd | Reel. 
k=2 kn? 


then 
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Since C includes B,* 


(1) lim P. = 0. 
For 2<k<£n, 
n—k+1 | 


1 
l=1 


Ax | 


lok 
n—k+1 


A; 
< : = | | 


n 
Ag dy 
l=k 


Therefore 


n n 
One| | < | | 
k=2 k=2 
= | | Vn—k+1 | | ue 
k=2 k=v+1 


where v may be chosen so that vy and »—v become infinite with 
n. Since lim, ..| V,/ | =0, for any e, y and m may be chosen suffi- 
ciently large so that 


(2) DY |] <e. 
k=2 
For 2<k<n-—1, we have 


n—2 
Vi | 


= 
C, 
l= dia | n—2 
2, 
< p+1 P | 


n—2 
On—141B1 
ist 


* Noérlund, Lunds Universitet, Arsskrift, (2), vol. 16 (1919), No. 3. 


= 
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where V,_, is the (n—2)th term of the sequence obtained when 
{| V, |} is summed by the matrix transformation t,», where 


n 
I=1 


for p<n, and t,,=0 for p>n. 

This transformation is regular, since A is regular; it follows 
that lim,...V,_2.=0, and that for any e, we may find n suff- 
ciently large so that 


tap 


(3) Ris é. 


From (1), (2), (3), and the fact that >" ,|u/| converges, it 
follows that lim,..|W,’|=0. This proves the theorem for 
v’=0. 

If V’0, we consider the sequence { V,— V’}; this sequence 


is summed by B to 0. Hence the Cauchy product of > wea by 
[v: — V’]+>—~_.o,. is summed by C to 0; that is, 


n 
>. 


lim| — = ®. 


Since C includes A, 


n 


lim k=1 
= 


therefore lim,...W,’ = U’V’. 


THeoreM 2. If is summable A to U', and if 
is summable B to V', then ).,_,w: is summable D to U'V'. 


PRroor. We have 


1 n n 
n k=1 k=1 


where for kn, and gn. =0 for k>n. 
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Since A and B are regular, this method of summation is regu- 
lar and limy..2n,n—%41=0; it follows that* 


k=l 
which completes the proof. 
For the proof of Theorem 3 we require the following lemma. 


Lemma. If {X,} and {[>°?_,Biyn—2+1]/Bn} converge to X and 
Y, respectively, and if { [>}_,beyn—z11]/bn} is bounded, then 
n—k+1 
( 
l=1 
{ = XY, 
k=1 
provided that (a) (b) <M for all n, 
where M is a positive constant; (c) T' includes T, where T’ and T 
are triangular matrix transformations defined by t and 
bnk 


PrRooF. Let 


n—k+1 
biVn—-k-142 
=y 
n—k+1, 

and let Let Z,= From (a), (b), and 
(c), it follows that lim, ..¢nt=0, <M’, where M’ isa 
positive constant, and lim,..)>.,-,¢nr= Y. Choose p such that 
for a given e>0, |X,—X| <e/(2M’) when k>p. For kSp, 
|X.—X| <L. Then 


n p n 
S coe ||Xe—X] + | Xe — X| 
k=l k=l k=ptl 


IIA 


P 
—— | |. 


Choose N > p, and such that | <¢/(2pL) for n> N. Then for 
n>N, |Z.—X <e, which proves the lemma. 


* Dale, American Journal of Mathematics, vol. 47 (1925), p. 82. 
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THEOREM 3. If }>*_,ux is summable A to U' and yi sum- 
mable B to V’, and if 
| bati +--+ + it, 


then >. ,ws is summable C to U’V’'. 


< 


Proor. Consider the triangular matrix definition 


A 


This definition satisfies the three conditions of the lemma, for 


A KOn—k41 A 


(4) < < 


Do 
l=1 l=k lak 


(5) >>| ena | = 1; 
k=1 
(6) A’ = CB’, 
where A’, B’, and C’ are triangular matrix definitions with 
Qnk = On,n—k+1 5 = ? Cnk = 


The definition C’ is regular. The theorem follows immediately 
from this lemma. 


THE GEORGE WASHINGTON UNIVERSITY 


n n 
n 
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ON THE LOCUS OF AN ANALYTIC EQUATION 
IN THE REAL PLANE* 


BY A. B. BROWN 


I have been urable to find in the literature a statement or 
proof of the following theorem. 


THEOREM. Let f(x, y) be analytic at the real point (xo, yo), with 
Sf (x0, vo) =0 and f(x, y) irreducible at (xo, yo).t Then the locus of 
the equation f(x, y)=0 in the real xy plane near (xo, yo) con- 
sists of one of the following three: (1) the point (xo, yo); (2) a single 
smooth curve through (xo, yo); (3) a cusp with vertex at (xo, yo). 


More detailed descriptions of (2) and (3) are contained in the 
proof which follows. 

By a change of coordinates we may suppose xp =yo=0. Ac- 
cording to the Weierstrass preparation theorem for the case of 
one independent variable, since f is irreducible at (0, 0), either 
f(x, y) =xQ(x, y), with Q(0, 0) #0 and Q analytic at (0, 0) or 
(1) f(x, y) = [y™ + +--+ + An(x)]Q(x, y), 
with 2 as above, m>0, A;(x) analytic at x=0, and A,(0)=0, 
(j=1,---, m). Since in the first case the real locus f(x, y) =0 
is merely a straight line, it is sufficient to consider the case that 
(1) holds. 

Since f is irreducible at (0, 0), the same is true of the algebroid 
function in (1), and hence its m-leaved Riemann surface is con- 
nected near (0, 0) and we can uniformize locally as follows: 

(2) x=", 
(3) =at+af+---, 


with y analytic at =0, and a neighborhood of the origin in the 


Il 


* Presented to the Society, February 23, 1935. 

+ That is, not the product, near (xo, yo) in the 4-space of the complex 
variables, of two functions each analytic and zero at (xo, yo). For theorems 
which we use involving functions of complex variables, see W. F. Osgood, 
Lehrbuch der Funktionentheorie, vol. 1, Chapter 8, §14, and vol. 2, part 1, Chap- 
ter 2, §§2, 4, 7,9, 10, 11. 
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t plane giving exactly the points (x, y) in a neighborhood of (0, 0) 
satisfying f(x, y) =0, each just once. Now x is real if and only if 


(4) gm (k = 0,1,---,m-—1), 


with 7 real and positive, negative, or zero. Substituting from (4) 
in (3), for a fixed value of k we get 


(5) y = = t+ beer? , = 


Either all the coefficients in (5) are real, in which case equation 
(4) for the given k, together with (2) and (3), gives only real points 
(x, y), or else there is a first non-real b;,;. In this case we have 
from (5) that 


y = R(t) + + be + +--+); 


where R(r) is real. Site b;,; is not real, we see that for suffi- 
ciently small 7 the parenthesis is not real. Since 7 is real, it follows 
that for small values of 70, y is not real. When 7=0, y=0. 
We thus see that in this case only r =0 gives real y. Next we show 
that at most one value of k can give real y for real 70. 

Suppose ¢t = and with 0<j—k<™m, both gave 
real loci for real 7, when used with (2) and (3). Then 
and b;,,=a,e7"*/™ would be real, (v=1, 2,---), and hence if 
a,~0, their ratio e-*)”*4/™ would be real. Therefore (j7—k)v/m 
would be an integer. Let m=m,-mz, where m, is the H. C. F. 
of (j—k) and m. Thus mz>1, since 0<j—k<m. Then v would 
have to be a multiple of mz, and hence 


(6) 


with ¢ analytic at the origin. Now if we set t:=¢e?*/™,then 
and ty"? =t™2e2*' =1™2, Hence from (2) and (6) 
we see that #4; and ¢ would give the same (x, y). But since m2 >1, 
t; ~t if t~0, and therefore we would have a contradiction. Con- 
sequently at most one of the m values of k gives real (x, y) when 
used in (2), (3), and (4) with real r 0. 

We now consider the case that some value of k does give a 
real locus. We define as the branches of the locus for the & in 
question the parts for which r=0 and 70, respectively. Thus 
the two branches have the origin as a common end point. From 
(4), (2), and (3) we see that on each of the two branches y 1s a 
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single-valued real function of x, analytic when x0. If m is even, 
x has the same sign on each branch. If m is odd, x has opposite 
signs on the two branches except at x =0, and y is a single-valued 
real continuous function of x for x on an interval of the real axis 
containing the origin. 

From (2) and (3) we see that (i) if a3=a.= --- =adn=0, then 
limz.oy/x=0; (ii) if ai=--- =dn1=0 but an~0, then 
lim,.o y/x =adm3£0; (iii) if some a;#0 with 0<j<m, then when 
x approaches zero, y/x becomes infinite. Hence the two branches 
have a common tangent line. By a rotation of axes we now ar- 
range that the x axis is the tangent line at the origin, and that 
x20 on the real locus if m is even. Since lim,.o y/x must then be 
0, we must have case (i), so that a;= --- =adm=0. Then, for 


7 dy = dy dt = m aks 1 
(7) = [(m + l)amyit™ + (——). 


Therefore lim... dy/dx =0. Hence we have a continuously turn- 
ing tangent; though if m is even, properly speaking we may only 
say that each branch has the positive x axis as right-hand tan- 
gent at the origin (a cusp). 

Now from (7), for +0, 


dy 
(8) — = bn+2t? 
dx 


Hence when ¢0, 


(0) d*y [-(2)] dt 1 ) 
dx? di \dx dx 


Therefore on either branch, with t#0, d?y/dx? is either iden- 
tically zero or nowhere zero. Hence unless the locus is a straight 
line, d*y/dx* has a fixed sign on each branch except at (xo, Yo) 
where it may be zero, and the branch does not meet T, the tangent 
at (xo, Yo), except at (xo, yo). We also see that as a point on the 
locus approaches (xo, yo), the curvature either approaches a fixed 
limit, possibly zero, or becomes infinite. The curvature at (xo, Yo) 
exists for each branch and equals this limit (also in the infinite 
case, where we have infinite curvature at (xo, yo)), as follows 
from (2), (8), and (9). 


= 
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As examples we mention y*—x'=0, y'—x'=0, and 
with the determination of (1—x?)1/? 
which equals 1 when x =0. These give respectively, at the origin, 
a minimum, a point of inflection, and a cusp with both branches 
concave upward. In none of the three cases is y analytic in x 
at the origin. An example where the locus is a single point is 
given by y+ix=0. 

In the case of a reducible function f(x, y), the real locus 
f(x, y)=0 neighboring (xo, yo) consists of a finite number of 
configurations of the kind described in the theorem, no two of 
which have any point except (xo, yo) in common. This is easily 
proved by use of theorems on resultants and on divisibility of 
one function by another. Of course two irreducible factors may 
have exactly the same locus. 


Co_uMBIA UNIVERSITY 


A PARTIAL DIFFERENTIAL EQUATION CONNECTED 
WITH THE FUNCTIONS OF THE 
PARABOLIC CYLINDER* 


BY HARRY BATEMAN 
The partial differential equation 


(av 
(1) ( 4+ =0, 
s==] 


0x2 OX, 


which was considered by Mehlerf in 1866, is a slight modifica- 
tion of an equation which occurs in wave-mechanics in the 
theory of the rotator in a plane and in space.{ The case in 
which » is a positive integer is then of chief physical interest 
and Mehler’s simple solution 


Pp Pp 
s=1 s=1 


acquires a physical significance. The function //,,(x) is the poly- 
nomial of Laplace and Hermite defined by the equation 


* Presented to the Society, December 2, 1933. 

+ F.G. Mehler, Journal fiir Mathematik, vol. 66 (1866), p. 161. 

t A. Sommerfeld, Atombau und Spektrallinien, wellenmechanischer Ergan- 
zungsband, 1929, p. 23. 


= 


1935-1 PARABOLIC CYLINDER FUNCTIONS 885 
(3) H,(x) = Dz)™e~* = (x), 


where 22 =x’. The H-notation is that adopted by Appell and de 

Fériet in their book,* the ¢-notation is one used in statistics. 
When p=2 we have the equation 


x—— y— + wW =0, 


(4) 
Ox oy 


which is closely related to one considered by E. L. Ince.t The 
normal or simple solutions are of type 


(5) V = hy_m(x)kn(y), 
where h(x), k(x) are solutions of the differential equation 
(6) 

dx? dx 


Ince generalized this solution by multiplying it by an arbitrary 
function f(z) and integrating with respect to 2. We shall use 
summation instead of integration and shall suppose m+1 to bea 
positive integer. We shall, moreover, put h,(x) =k,(x) =H,(x), 
where [H,(x) is defined for one-half of the x plane and for all 
values of s by the equation 


‘| 
-{ + + (x — y)]dy, 


which reduces to the well known formulaof Laplacet and Mehler§ 
when m=0 and s is a positive integer. On the right-hand side 
of the equation (7), the many-valued function (x + zy)" is sup- 
posed to be defined for y>0 so that when x—0 it reduces to 
y?"exp(+in7r). When x is complex and unrestricted the inte- 
gral is not single-valued, but if x is restricted to the half-plane 


* P. Appell and J. Kampé de Fériet, Fonctions Hypergéométriques et Hyper- 
sphériques, Polynomes d’ Hermite, 1926, p. 344. 

7 E.L. Ince, Proceedings of the Royal Society of Edinburgh, vol. 44 (1924), 
p. 242. 

P.S. Laplace, Théorie A nalytiques des Probabilités, Livre 2, Oeuvres, vol. 
7, p. 299. 

§ F. G. Mehler, loc. cit. 


— 
= 
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R(x) >0, the integral agrees in value with a single-valued func- 
tion which satisfies the differential equation (6). It is, indeed, 
easy to verify by an integration by parts that the definite in- 
tegral has the form on the left, where H,(x) is some function of 
s and x such that H,(x) satisfies (6). The proof of the last prop- 
erty depends on the fact that (x+7y)’ is a solution of the differ- 
ential equation (4). 

To find which solution of (6) the integral represents we first 
choose R(s-+-m) >0 and make x-0. Writing s = 27 and observing 
that* 


(8) f 


we find that 


x 


1 1 


H.,(0)r : = 
(> (;): 


(0)T(— r) = (— =). 


(9) 


Hence, with the usual notation for the confluent hypergeomet- 
ric function (without the suffixes), 


2 2°2 
1 


(10) 


In Whittaker’s notation, if s is not a negative integer, 


(11) H,(x) = exp (=) Du. 


This equation may be used to define the function H,(x) for all 
real and complex values of x. An equation equivalent to this has 


* The value of the integral is given in a slightly different form by G. N. 
Watson, Proceedings of the London Mathematical Society, vol. 8 (1910), p. 
393. 


— — — 

2 2 
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indeed been used for this purpose* for all real and complex 
values of s except negative integral values. When R(s-+m) >0 
we can, by integrating the integral in (7) a sufficient number of 
times, pass to the case in which R(s+m) <0. An exceptional 
case arises when s is a negative integer. If we take m+s to bea 
negative integer, the integral has a definite value and gives a 
definition of H,(x), but a different notation is used for the func- 
tion thus defined.{ If m-+s is a positive integer, the integral and 
Cn+s,m are both zero for s<0. When s is a positive integer 
p—n, the relation (7) gives at once the expansion 


Pp 
(12) iy)? = Cy.mi™H m(y), 

m=0 
which may be checked by comparing the terms of degree p on 
the two sides of the equation and noting that both sides satisfy 
the partial differential equation (4). This equation may be used 
to derive from (7) the expansion (s not a negative integer) 


Pp 


(13) = (— 1)" p-m(*), 


m=0 


which may be inverted so as to give an expansion for the product 
in the formt 


(14) H,(x)H,(x) = 


m=0 


This relation may be established directly by induction.§ When 
s =n the equation takes the form 


(15) [H,(x) |? = 
s=0 


* M. Plancherel and W. Rotach, Commentarii Mathematici Helvetici, vol. 
1 (1929), p. 227. See also S. C. van Veen, Amsterdam Proceedings, vol. 34 
(1931), p. 257. 

+ G. N. Watson, loc. cit.;O. Volk, Mathematische Annalen, vol. 86 (1922), 
p. 296; Appell and de Fériet, loc. cit., p. 362. 

t An equivalent expansion is given by Gorakh Prasad for positive integral 
values of s and m. Products of Hermite polynomials are considered also by 
N. Nielsen, Recherches sur les Polynomes de Hermite, Copenhagen, 1918. 

§ Gorakh Prasad, Proceedings of the Benares Society, vol. 2 (1920), p. 18. 
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This series may be compared with Mitra’s series* for D,? (x), 

(16) [D,(x) ]?= 
s=0 


which has been written in an abbreviated form with the aid of 
Schuster’s notation 


(2s)!! = (2s)(2s — 2)--- (2), 
(2s + 1)!! = (2s + 1)(2s — 1)--- (4), 
and a natural extension of it, namely, 
(— 1 — 2s)!"(— 1)(— 3)--- (1 — 2s) = 1. 


Mitra’s expansion may be generalized with the aid of the well 
known recurrence relations 


(17) Da+i(x) Tt xD,(x) + nD,,-1(%) = 0, 
1 

(18) Di (x) + xD, (x) — nD,_1(x) = 0, 

(19) Hayi(x) — xH,(x) + nH,_1(x) = 0, 

(20) A! (xs) = 


The resulting equations are 
(21) “ 
= (2r + 2s — 1)!!(2n — 2r — 2s — 1)!"(2s — 2r — 1)!A,, 
s=0 
where 
A,-(2s)!(2s — = Ae,(x); 
(22) 
= (2r + 2s — 1)!(2n — 2r — 2s — 1)!(2s — 2r + 


s=0 


where 


* S.C. Mitra, Proceedings of the Edinburgh Society, (2), vol. 4 (1934), p. 
27. 

+ A. Schuster, Philosophical Transactions, Royal Society of London, (A), 
vol. 200 (1903), p. 181. The extension is also given by Schuster. 
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C,(2s)!(2s 4+ = 


These equations may be checked* by comparing them with (13) 
and (14) after these have been multiplied throughout by D¢ (x) 
so that they become equations between D-functions. The verifi- 
cation depends on the identities 


(2m + 2n — 1)!!(— 2m — 1)!(2r — 1)! 
= 2*Cm.eC minsr.eS! (2m + 2r — 2s — 1)! 


s=0 
(23) - (2s — 2m — 2n — 1)!!(2m + 2n — 2s — 1)!!; 
(2m + 2s — 1)!!(2m — 2s — 1)!1(2s — 2m — 1)! 


= D> (— 2)"Cp,mCon—p,mm!(2p — 2n + 2s — 1)!! 


m=0 


- (2n — 2m — 2s — 1)!!(2s + 2n — 2p — 1)!!. 


Let us now consider the integral 
(24) Hy,m(a, x) = f cos a+ y sin a)o(™(y)dy, 


which is convergent when —17<2a<z if H,(z) is defined by 
(11), for the behavior of D,(z) when z is large is known. Since 


— H,(z) = vH,_1(2), 
02 

we readily find that 


x) = — x sin vH,_1( )o(y)dy 


(25) + cos. (yay 
= — yx sin aH,_1,0(a, x) 


+ v(v — 1) cos @ sin aH,_2,0(a, x). 


* The expansions obtained by putting n=2r in (21) and n=2r—1 in (22) 
have already been given by R. K. Varma in a slightly different form, Proceed- 
ings of the Benares Society, vol. 10 (1928), p. 11. Some special cases were con- 
sidered previously by Gorakh Prasad, Proceedings of the Benares Society, 
vol. 2 (1920), p. 12. 


| 
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Since, when a—0, 


(26) H,,o(0, x) = H,(x) = xH,_1(x) — (v — 1)H,_2(x), 
we find that a solution of equation (25) is 
(27) H, o(a, x) = cos’aH,(x). 


Integrating the integral (24) by parts m times we find 
(28) Hi, (a, %) = C,,_m! cos’—™ a sin” a H,_»,(x). 
The associated equation* 
H,(«% cos a + y sin a) 


(29) 
m=0 


was obtained by Kampé de Férietf in the case when 1 is a posi- 
tive integer by using the generating function of the Hermite 
polynomial, a method which he used also to obtain an expansion 
for H(pixitpexe+ -- - +,x,) in a series of simple functions of 
type (1) on the supposition that p?+ 2+ --- =1. In the gen- 
eral case when 7 is not a positive integer we may consider the 
series in powers of tan a of sec’ aH,(x cosa+ysina). Ifvisnota 
negative integer, the radius of convergence of this series is de- 
termined by the location of the singularities of sec’ a, for H,(z) 
is uniform and analytic over the whole of the finite part of the 
z-plane. The power series for sec’ aH,( ) is thus convergent for 
—a1 <2a<7, and so to establish the truth of the series under this 
condition it is only necessary to prove that if ¢=tan a, then 


I(m, v) = |— sec’ aH,(x cosa + y sin a) | 
(30) dt™ t 


ll 


= C, 


* When » is a positive integer this equation may be checked by observing 
that both sides are polynomial solutions of (4) with the same terms of degree 
vy and these terms determine the solution uniquely. 

+ J. Kampé de Fériet, Danske Videnskabernes Selskabs, Mathematisk- 
Fysiske Meddelelser, vol. 5 (1923), No. 2. See ai.. Appell and de Fériet, loc. 
cit., and V L. Mutatker, Journal of the Indian Mathematical Society, vol. 1 
(1934), p. 53. The case in which a=7/4 is due to C. Runge, Mathematische 
Annalen, vol. 75 (1914), p. 130. The equation was in fact, used by him to solve 
a special type of non-linear integral equation. 


| 

| 

| 

| 
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This proof may be carried out by induction. Thus, if we suppose 


that sec’ aH, (x cosa+ysin a) = F(v,‘t), we have, on differenti- 
ating first once and then m times, 


(y — xt)yF(v — 1, t) + viF(y, 


d 
(1+ #) — Fs, 
dt 


I(m + 1, v) = yol(m, v — 1) — maxvI(m — 1, vy — 1) 


+ m(v — m — 1)I(m — 1, »), 


for ¢=0, and the recurrence relation is satisfied by the expression 
on the extreme right of (30). 

When vis a positive integer, the integral representing H, (a, x) 
may be evaluated by putting z=xcosa+vysina and using 
Mehler’s expansion* 

csc a-exp [esc? a(x? cos? a — 2x2 cos a + 2? cos? a) | 
(31) | 


s=0 $3 


When the integral for H,,»(a, x) is treated in the same way, we 
obtain the integral 


(32) f H,(2)om(z csc a — x ctn a)dz 


= C,, mm! cos’—™ a aH,_m(x), 
which gives the expansion 


csc a@ — x ctn a) 


(33) = (24)—1/2 > — m)!]-! cos*-™ sin™ aH,(z)Hn—m(x). 


Another expansion may be obtained by writing ctn a=e~* and 
observing that the function V =csc’ a-H,(x cos a+y sin @) isa 
solution of the partial differential equationT 


* F. G. Mehler, loc. cit. 

+ This equation may be transformed into one which Laplace solved by 
means of a definite integral, Théorie Analytique des Probabilités, Livre 2, 
Oeuvres, vol. 7, p. 294. The equation was used later by Smoluchowski in his 
theory of the Brownian movement. See B. Hostinsky, Annales de 1I’Institut 
Poincaré, vol. 3 (1932), pp. 1-72. 


| 
| 
| 
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(34) — = — 
Ox? Ox at 


If this equation has a solution of the form 
(35) V = Cilt)x*/k!, 
k=0 


the coefficient C;,(¢) must satisfy the differential-difference equa- 
tion 
(36) Cé (t) + = 
and if the value of C,(0) is known, an appropriate solution is 
(37) Cr(t) = — 

s=0 


Instead of defining C;(¢) as a coefficient we may define it by the 
equation 


a'V 
to = 


In particular, when V =csc’ a-H,(x cos a+y sin a), we have 
(39) C.(t) = C,,4 k! cos* a csc” aH,_x(y sin a), 
and if n=v—k, the foregoing expansion for C;(t) leads to the 
equation 

H,(y sin a) 


40 


s=0 
which may be rewritten in the form 
(41) H, (uv) = >> — 
s=0 


When 2 is a positive integer, this expansion is equivalent to one 
given by Appell and de Fériet.* The associated integral 


* Appell and de Fériet, loc. cit., p. 346. Equation (40) is related also to 
equation (67) in the paper of S. Goldstein, Proceedings of the London Mathe- 
matical Society, (2), vol. 34 (1932), p. 103. 
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(42) f H = m!H,(0)v™—-*(1 — v?)*/2/s! 


is closely connected with one given by Doetsch.* 

If in the expansion (29) we put a=7/4 and set first x = u(2)'/?, 
y =v(2)*/2, secondly x =(u+v)(2)!/2, y=0, we obtain two ex- 
pansions for H,(u+v). Equating them and writing x for u(2)*/?, 
y for v(2)'/?, we obtain the equation (for vy >0) 


m=0 m=0 


This was obtained by Varmaf for the case in which 1 is a posi- 
tive integer. When x =y Runge’s expansion 


2°H,(x(2)!/2) = Cy mH m(x) 
m=0 


suggests that it may be profitable to study the expansion of a 
product H,(x)H,(x) in a series of functions of type H,(x(2)*/?). 
Such expansions have been studied by Mitra.f{ 


CALIFORNIA INSTITUTE OF TECHNOLOGY 


* G. Doetsch, Mathematische Zeitschrift, vol. 32 (1930), p. 587. 

+ R. S. Varma, Proceedings of the Benares Society, vol. 9 (1927), p. 31. 

S.C. Mitra, Bulletin of the Calcutta Mathematical Society, June, 1926; 
Proceedings of the Benares Mathematical Society, vol. 9 (1927), p. 21; Pro- 
ceedings of the Edinburgh Mathematical Society, (2), vol. 4 (1934), p. 27. 
See also S. C. Dhar, Journal of the Indian Mathematical Society, new ser., 
vol. 1 (1934), p. 105. 
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A GENERALIZATION OF THE BERNOULLI 
POLYNOMIAL OF ORDER ONE* 


BY B. F. KIMBALL 


1. Introduction. The idea of a Bernoulli polynomial, general- 
ized so as to allow of a continuous variation of the index of its 
degree, is not a new one.f In the present paper the writer has de- 
veloped a simple definition of a generalized Bernoulli polynomial 
of order one which brings the Bernoulli polynomial into direct 
relationship with the generalized Riemann zeta-function. An- 
other very interesting property of this generalization of the 
Bernoulli polynomial is brought out in §7. 


2. Difference Equation Considered. For the purposes of this 
paper in dealing with a complex power s of a complex number 
t, if t=re®, then =erllog <9 <7). The following 
difference equation will frequently be referred to: 


1 
(1) i [f(x + w, s) — f(x, s)] = sx, 


where the difference interval w is taken real and positive. Also 
there will be occasion to impose the asymptotic condition: for any 
value of s such that R(s) <0, 


(2) f(%,s) as R(x)> +0. 


There will also be occasion to refer to the following regions on 
the x and s planes: 

Region X. All points on the x plane other than the negative 
axis of reals and the origin. 

Region S. All points on the s plane within and on the bound- 
ary of a circle of radius M centering at the origin, (usually 
taken arbitrarily large, see §5). 


3. Uniqueness of an Analytic Solution. 


* Presented to the Society, June 23, 1933. 
t See N. E. Nérlund, Vorlesungen iiber Differenzenrechnung, p. 53. 
Future references to this book will be indicated by the letter N. 
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LEMMA 1. Given a function of x and s which is defined for all 
values of x on region X and which for any such value of x is analytic 
in s over region S. Given also that this function is a solution of the 
difference equation (1) for x in X and that it satisfies the asymp- 
totic condition (2). Such a function is unique (the existence of 
such a function is established in §5). 


Proor. Assuming the existence of two such functions f1,2(x, s), 
we have 


(3) w,s)—¥(x%,s)=0, (xeX;~=fi — fr). 


From (2) it follows that ¥(x, s) must be identically zero, xc X, 
if R(s) <0. Let xo denote any value of x in region X. We have 
then from the hypothesis of Lemma 1 that ¥(xo, s) is an analytic 
function of s over the region S. Thus this analytic function must 
be zero for all values of s in S, and the lemma follows. 


4. Definition of a Generalized Bernoulli Polynomial. For s=a 
positive integer or zero and x real, the classical Bernoulli poly- 
nomial, B,(x) =B,(x, 1), has been defined (N., p. 18) as the poly- 
nomial solution of the difference equation (3), with difference 
interval equal to unity, which takes on the value of the corre- 
sponding Bernoulli number B, when x =0. Nérlund has termed 
the polynomial above the “Bernoulli polynomial” (of order one), 
with difference interval equal to one, and has defined Bernoulli 
polynomials of positive and negative order m for all integral val- 
ues of m, where the difference interval w is not restricted to be 
equal to unity (N., pp. 129, 138). 

The author has found it of interest to set up a definition of 
a generalized Bernoulli polynomial B,(x, w) of order one where 
s and x are allowed to take on complex values. To avoid compli- 
cations, w has been taken as real and positive. 


DEFINITION. The solution of the difference equation (1) de- 
scribed in the statement of Lemma 1 shall be considered as the 
generalized Bernoulli polynomial B,(x, w) of order one for s in 
region S. 


5. Solution Based on Euler-Maclaurin Summation Formula. A 
solution based upon the Euler-Maclaurin summation formula 
will first be set up. We define a periodic function B,,(#) as fol- 
lows (N., p. 30). Let the variable ¢ be restricted to real values 


= 
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and denote by 3,,(#) the classical Bernoulli polynomial of de- 
gree m. The function B,,(é) is taken so that 


(4) B,.(t) = Bn(t), (0<t <1), 


and B,,(t+1) —B,,(t) =0 for all values of t. Let m be any positive 
integer greater than a positive number M (radius of circular 
region S). Define the function f(x, s)=F(x+hw, s) by the rela- 
tion (see N., p. 52): 


F(x + hw, s) = 
(5) k=0 


+ (m + f — z)(x + wz)*-™"dz, 
<h<1). 


It is not difficult to verify the following statements concerning 
this function for x in region X and s in region S. 


(i) For any fixed value of x in region X, f(x, s) exists and 1s an 
analytic function of s in S. 

(ii) The function f(x, s) defined above ts a solution of the differ- 
ence equation (1). 


This is easily verified as follows. Consider the infinite in- 
tegral as a sum of an infinite series of integrals with limits at 
0, w, 2w, - - - . Then take the difference of each side of the equa- 
tion with respect to x with difference interval w. This reduces 
the infinite series of integrals to one integral, and the expression 
on the right is easily recognized to be the Euler-Maclaurin ex- 
pansion of s(x-+/w)*—! in terms of its first differences: 


(iii) For any value of s such that R(s)<0, f(x, s)—0 as 
R(x) 

(iv) The function f(x, s) is the generalized Bernoulli polyno- 
mial, B,(x, w), of order one for s in region S. 

(v) If w=1 and s ts a positive integer (or zero) less than M, 
f(x, s) reduces to the classical Bernoulli polynomial. 


(vi) The above region of definition S of B,(x, w) on the 
s plane may be made as large as desired by giving to m a suffi- 
ciently large value in equation (5). Thus we have here a proof 


= 
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of the existence of B,(x, w) and its analyticity with respect to s 
over the finite s plane forxc X. 

(vii) Set h=0 and x=w in (5). This brings out the relation 
(N., p. 129) 
(6) Bw, w) = w*B,(1, 1) = w*B,(1). 

6. Multiplication Equation Satisfied by* B,(x, w). The differ- 
ence equation 

h(x + w/n) — h,(x) = wsx*, (x¢X), 


is satisfied by Bx +kw/n, w) and n'*B,(nx, w). When 
R(s) <0 the asymptotic value of each of these functions, as R(x) 
becomes positively infinite, is zero. Thus, with the reasoning 
employed in the proof of Lemma 1, it follows that 


n—1 


(7) n'-*B,(nx, w) = B(x + kw/n,w), (xeX), 


k=0 


for all values{ of s. In particular for n=2 and x =w/2, this gives 
the relation 


B,(w/2, w) = — 1)B,(w, w), 
which taken in connection with (6), leads to 
(8) B(w/2, w) = w*B,(1/2). 


7. An Interesting Property of B,(x, w). Define function ¢,(x) by 
the relation: 


(9) ¢.(x) = B,(w/2 — x, w) — eT*B,(w/2 + x, w). 


Using the fact that B,(x, w) is a solution of the difference equa- 
tion (3), one arrives at the conclusion that 


(10) $.(% + w) — ¢,(x) = 0, (xe X). 


When s is a positive integer or zero, it follows from well known 
properties of the classical Bernoulli polynomials that ¢,(x) is 
identically zero for all x. Thus, (i) ¢,(x) is a periodic function 


* See Thiruvenkatacharya, On some properties of the zeta function, Journal 
of the Indian Mathematical Society, vol. 19 (1931), pp. 92-96. 

Tt In §5 the existence of B,(x, w) and its analyticity with respect to s over 
the finite s plane was proved. 
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with period w for all values of s; (ii) when s is a positive integer 
or zero, @,(x) is zero for all values of x; (iii) ¢,(0) is given by 


$.(0) = (1 — e**)B,(w/2, w) = (1 — e***)w*B, (1/2). 


These properties point to a significant role played by B,(1/2) in 
the behavior* of w). 


8. A Closely Related Integral. Consider the integral} L(x, s). 


w 2 
(11) = («+ dz, 
Qrid _in w cos 
(R(x) > 0). 


We may write L(x+w, s) in the form 


w w+ ico 2 
w-iw w cos 


The integrand of these integrals is analytic as a function of z 
over the region bounded by the two lines R(z) =0 and R(z)=w 
except for a pole at z=w/2. The residue of this pole is 
s(x-+w/2)*-1. Hence the integral taken around the rectangle 
bounded by the two lines referred to above, and lines parallel 
to the axis of reals at —iN and +iN (N real and positive) comes 
to ws(x+w/2)*. As I(z)—>+ «, the absolute value of the inte- 
grand goes to zero like f*e~*', i>», uniformly for 0< R(z) <w 
(c real and positive). It follows that if we let N— «, the contri- 
butions to the circuit integral made over the above lines at 
+iN approach zero. Thus, from (11) and (12), we have 


(13) L(x+w,s) — L(x, s) = ws(x + w/2)*", (R(x) > 0). 


Moreover, this relation is valid for all complex values of s. Also 
it is not difficult to show that for R(s)<0, L(x, s)--0 as 
R(x)—> «. Upon applying the reasoning used in proving Lemma 
1, it follows that L(x, s) is the unique solution of difference 
equation (13) for x on the right half-plane which satisfies the 
asymptotic condition (2). Hence, for x on the right half-plane, 
L(x, s) is identical with the generalized Bernoulli polynomial 
B,(w/2+x, w). 

* See N., p. 130. 

T See N., p. 75. 
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9. Relation to the Generalized Riemann Zeta-Function. The 
generalized Riemann zeta-function may be introduced by the 
relation 
(14) f(s, 1 s) 

x,1—s)= 
(% + 
For R(x)>0, the definition of ¢(x, 1—s) as a function of s has 
been extended to include the finite s plane and under such defini- 
tion s{(x, 1—s) is found to be an entire function of s.* 

On the other hand, from the difference equation satisfied by 
B,(x)=B,(x, 1), it follows that 


(xe X, R(s) < 0). 


(15) > : [B 1) — B,(x)] 
If R(s) <0, it is known that B,(x+r+1)—0 as ro. Hence, 
1 
16 «eX, R(s) < 0), 
(16) EX RW) <) 
(17) Bx) = — sf(x,1—s), (xe X, R(s) < 0). 


When R(x)>0, the functions on each side of this relation, con- 
sidered as functions of s, are defined and analytic over the finite 
s plane. Thus the above equation holds for all values of s, and 
we can state the following theorem. 


THEOREM. For R(x) >0 and for all values of s, the generalized 
Bernoulli polynomial of order one is related to the generalized 
Riemann zeta-function by equation (17). 


As far as the writer is aware, this is the first time that the 
relation (17) has been set up for values of s other than integers. 
The fact that it has been worth while to study the zeta-function 
as a function of a complex variable, may, by reason of equation 
(17), mean that the study of a generalized Bernoulli polynomial 
will yield further interesting and worth while results. In passing 
it is interesting to note that the relation (17) may serve to ex- 
tend the domain of definition of ¢(x, s), for all values of s, to val- 
ues of x over all of region X. 


* See, for example, Wittaker and Watson, Modern Analysis, 3rd ed., 
Chap. 13. 
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From the relation between the zeta-function and the general- 
ized Bernoulli polynomial there emerge several interesting con- 
sequences. From the well known functional relation for ¢(1, s) 
it follows that 


(18) sBy_,(1) = sin (sx/2)T(2 — s)B,(1) 

for all values of s. Since B,(1/2) plays a central role in the gen- 
eral behavior of the function B,(x, w), it is of interest to apply 
the above relation to B,(1/2). From the multiplication equation 
(6) when we sct w=1, n=2, x=1/2, we obtain 

(19) B,(1/2) 1)B,(1), 

which, when substituted in (18), gives 

2*—1 
1 


(20) By,(1/2) = | |re — s) sin (sx/2)B,(1/2). 
This relation holds for all values of s. 

If w were given other values than unity, a further extension 
of the zeta-function would result from the relation (17). This 
would mean that in the series definition of the zeta-function 
(14) the denominators x+k would be replaced by numbers 
x-+kw. Such a generalized zeta-function would be the solution 
of the difference equation 


1 
(21) (eX), 
w 


analytic in s over the finite s plane (except for a pole at s =1) 
which for R(s)>1 has its asymptotic value equal to zero for 
R(x) positively infinite. 


AvBany, N. Y. 
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THE METHOD OF MOMENT DISTRIBUTION FOR 
THE ANALYSIS OF CONTINUOUS STRUCTURES 


BY LOUIS BRAND 


1. Introduction. If ab is any member of constant cross-section 
forming part of a continuous structure, the moments at its ends, 
Ma and Mia, are given in terms of the angle-changes 0,, & at 
its ends and the lateral deflection per unit of length R by the 
slope-deflection equations* 


Ma = 2EK + 3R) + 
Moa = 2EKy.(26, 3R) Ces 


Here K.s = Ky. denotes the sectional moment of inertia of ab di- 
vided by its length (//1) and Cu, Coa are the numerical values 
of the fixed-end moments due to the loading on ab. For the 
derivation of these equations and the sign conventions em- 
ployed, reference may be made to the Bulletin just cited. 

When there is no lateral deflection, or when this is neglected, 
R=0, and we write the slope-deflection equations 


(1) Mas = + 6s) + Mar, 
(0) 


(2) Moa = + 92) + Mia , 


in which M3), M@ denote the fixed-end moments inclusive of 
sign. Suppose also that at all joints of the structure, other than 
certain fixed ends, there is no external momental load. At any 
such joint a, we must have 


(3) > Ma = 0, 


the summation ranging over all members that meet at a. At a 
fixed end c, 8, is given. Owing to the continuity of the structure, 
all members meeting at a joint a rotate through the same angle 
6,. At each such joint we have an equation of type (3). Thus we 
have precisely as many equations (3) as we have unknown 
angles, so that in general these equations determine the angles 


* Bulletin No. 108, Engineering Experiment Station, University of Illinois, 
p. 20. 


= 
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uniquely. Having found them, equations (1) and (2) give the 
terminal moments for each member of the structure. 

Hardy Cross has devised a method of successive approxima- 
tions for solving equations (1), (2), and (3) for the terminal 
moments M,, that is very simple and effective. In his original 
paper* he regards the method as the mathematical parallel of a 
physical process. To quote Cross: “The beams are loaded or 
otherwise distorted while the joints are held against rotation; 
one joint is then allowed to rotate with the accompanying dis- 
tribution of the unbalanced moment at that joint while the re- 
sulting moments are carried over to the adjacent joints; then 
another joint is allowed to rotate while the others are held 
against rotation; and the process is repeated until the joints are 
‘eased down’ into equilibrium.” 


2. The Successive Approximations. To begin the calculation 
corresponding to this process we write at the ends of each mem- 
ber the fixed-end moments corresponding to the load it carries. 
Consider a joint a at which the members ab, ac, - - - , ak meet. 
If the algebraic sum of the fixed-end moments at a, 

is not zero, we add a moment —S at a and distribute it among 
the members meeting there in the proportion of their stiffness. 
The stiffness of any member abd is defined as the moment which 
must be applied at its end a to make 0,=1 when 3 is held fixed 
(@,=0). From (1) it is clear that the stiffness of ab is 4EK«s. 
Thus the stiffnesses of the members ab, ac,---, ak are pro- 
portional to Kas, Kac, --- , Kaz. The balancing moment —S is 
therefore distributed among these members so that they receive 
respectively 

Kas Kat 

-S§,- = -§ 
>> Kai > Kai 

While joint a is allowed to rotate all other joints are held 

fixed: in particular, #,=0.= - -- =@,=0. Hence from equations 


* Proceedings of the American Society of Civil Engineers, May, 1930, pp. 
919-928. For the paper with the complete discussions, see the Transactions of 
the American Society of Civil Engineers, vol. 96 (1932), pp. 1-156. 
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(1) and (2), a moment 4£K,,0, applied at end a of ab induces a 
moment 2EK;,,@, at the fixed end b of ab. Thus, in the case of 
members of constant section, one-half of the moment distrib- 
uted to ab in balancing the joint a is carried over to the end b. 
We therefore call 1/2 the carry-over factor for members of con- 
stant section. 

The three basic operations involved in the above method of 
moment distribution are therefore: 

BD (Balance and Distribute). If the moments at a joint do 
not balance, add the balancing moment at the joint and distribute 
it to the members meeting there in the proportion of their stiffness. 

C (Carry-over). One-half of each moment distributed to a mem- 
ber is carried over to its other end. 

The moments carried over in operation C destroy the balance 
achieved in operation BD. Hence BD must be repeated on the 
moments carried over. A succession of operations C and BD 
will be called a cycle. If, at the end of any cycle, the fixed-end 
moment at end a of member ab is added algebraically to all the 
distributed and carried-over moments there, the sum will give 
an approximate value of Mw». We shall show that this approxi- 
mate value approaches the true value of M. as the number of 
cycles performed increases indefinitely. 

We next consider in detail this process of computing 1... The 
first cycle consists in writing the fixed-end moments at a and 
distributing the balancing moment at a according to operation 
BD. As regards the end a of ab, this gives the entries 


MY? 
(1) Kav (0) 
BD1 Ma = — Mai . 
(BD1) 


The moment M{, distributed to ab at 6 must now be carried 
over to a according to operation C. All these moments carried 
over to a create a new unbalance at this joint, which is then 
balanced and distributed according to operation BD. Thus the 
next two entries for the end a of ab are 


1 
(C1) = Mie, 
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1 Ka 
(BD2) = where My : 
Similarly the next cycle of operations gives the entries at end 
a of ab: 


(C2) 


1 
(BD3) Ma, where « 


Continuing the cycles indefinitely and adding all the entries, 
we obtain 


Ma = us +[ a +— 
(4) 


(2) (3) 


If the end a of ab is fixed, equation (3) does not apply and 
hence the operation BD is omitted. In this case M is built up 
from M3 and the successive moments carried over from b, and 
the first series in (4) is absent. 


3. Convergence of Series. We shall first establish the conver- 
gence of the series in (4). Let C denote the sum of the absolute 
values of all the fixed-end moments for the entire structure. Then 
writing pa =Ka/> Kai, we have from (BD1) 


| Mes’ | paw | Mer | <C. 


Now at any joint a, other than a fixed end, 


MP | = (part poet + | Mee 


Forming such inequalities for all joints which are not fixed and 
adding the results we obtain 


> | My |< > | M 
Hence from (BD2) 


el < 


| Mo | < po | M 


— 
— 
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In the same way we may show that 


(n) 


|Ma|<cC, (n = 3,4,---). 


Both series in (4) are therefore absolutely convergent. 


4. The Series Satisfy the Slope-Deflection Equations. If we 
identify the first series in (4) with 4EK.@., then the second 
series must represent 2E Kia, and (4) reduces to equation (1). 
That this identification is justified is readily seen from the fact 
that after any operation BD, the summed moments at a joint @ 
(which is not fixed) satisfy equation (3); for the operation BD 
was expressly designed to accomplish this. Hence in the limit 
the moments M,, given by (4) satisfy the system of equations 
of type (3)—one equation for each free joint. Hence we see that 
the sums of the convergent series 


Ms; +— +---, (inota fixed joint), 


satisfy precisely the same system of linear equations that 
4EK;,0; satisfy. Thus if the system of equations (3) has a unique 
solution for the angles 0;, this solution is given by 


5. Moment at a Fixed End. If a is the fixed end of a member ab, 
the equation (3) does not apply and the operation BD is there- 
fore omitted. The fixed-end moment MQ and the moments car- 
ried over from 6 will total up to the actual external moment act- 
ing on the structure at a. For on putting 0.=0 in (1) we have 

(0) 
(6) Mw = 2EK a + Mar ; 
and from (4), on omitting the first series due to the operation 
BD and retaining the second due to operation C, 


1 1 


In view of (5) the series in brackets sums to 4EKi.4%, so that 


906 LOUIS BRAND [December, 


the limiting value of the right-hand side has precisely the value 
given in (6). 


6. Procedure for a Member Hinged at One End. If a is the 
hinged end of a member ab, the end of a cycle will always reduce 
the total moment at a to zero. But this series of operations BD 
and C leading to the result zero may be avoided by taking the 
stiffness of ab proportional to Ki =(3/4)Kas, and replacing the 
fixed-end moments M2, Mf at the outset by 


1 
=0, = > Mav. 


These are the terminal moments for ab due to the given loading 
when a is hinged and b fixed. No moments need then be carried 
over to a from b, and consequently there is nothing to carry back 
from a to b. To justify this procedure, we have from (1) and (2) 
in this case, 
0 = + 0) + Ms, 

Mie = 2EKax(20s + 02) + Moe 

On eliminating 0, from these equations we obtain 


3 1 
= 2EKas + Ms 


= Hy. 


But from the successive approximations for obtaining My,, with 
H{ as a starting moment and no carry-overs from a, we have 


1 1 

(0 (1) (2) (3) 

My. = + + > Moa + Ma. | 
where, in view of (5), the series in brackets is precisely 4EKi.9 
as all distributions to ba are now made proportional to Kia. 
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ON CERTAIN HIGHER CONGRUENCES* 
BY LEONARD CARLITZ 


1. Introduction. This note is concerned with the higher con- 
gruence 


(1) II ¢-@=A4 (mod P). 


deg G<m 

Here A, P, G denote polynomials in an indeterminate x with 
coefficients in a Galois field GF(p") of order p*. The product 
in the left member extends over all G of degree less than some 
fixed m; the modulus P is assumed irreducible of degree k. As 
will appear below, we may without loss assume k>m. 

The congruence (1) has either no solution at all, or else has 
p”” distinct solutions; if ¢ is any solution, then the general solu- 


tion is furnished by t+G, where deg G<m. Define o; by means 
of 


(4 + x)(u + x) + 


Put 


ll 


t+--- +e, 
P! = + (Rk — +--+ + 
Then we prove the criterion: The congruence (1) is solvableif and 


only if each product (o;/(Fm—1)”")AP’, (j7=0,--- , m—1), is con- 
gruent (mod P) to a polynomial of degree <k—1. 


2. Some Properties of ~m(t). We denote by Wm(t) the product 
appearing in the left member of (1). Also, we let 


m—1 m—1 

Fn = 2"), Ln = — x), Fo= =1. 
im im 

Thent 


* Presented to the Society, September 10, 1935. 
+ L. Carlitz, Duke Mathematical Journal, vol. 1 (1935), p. 141. 
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ni 
i=0 


(2) van) = 92 ) 


We shall require a generalized form of identity (2). From the 
definition of y,,(Z), it is clear that 


(3) = II + + Cmit™ + --- + 69), 


the product extending over all sets (cn, - - - , co) in the GF(p"). 
The right member of (3) may be written in the form 


TI IT’ + + + 
+ + co) } 


the inner product extending over (Cm, , €j41, €j-1, °° , €o) 
only. If we put 


= TT + + + cp +--+ +09), 


then (3) becomes 


(4) Vmii(t) = II fit + cx?). 


Now by a formula due to E. H. Moore,* the polynomial f;(¢) 
may be written as a quotient of two determinants: 


A ;(t) 
(5) fit) = 
D; 
wheref 
| 
where i=0, 1,---, m, and 


* This Bulletin, vol. 2 (1896), p. 191. 
+ For brevity we write only the elements in the ith row: | aio aiy-** ais| : 
(i=0, i, S). 


ce 
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where 1=0, 1, - - - , m—1. By (5) and (6) it is clear that 
Silt + cx?) = fi(t) + of 
We shall evaluate f;(x’). Using (6), we have 


(7) (i=0,---,m), 


The calculation of D; is somewhat more involved. Expanding 
the determinant 


with respect to the elements of the first row, we have on the one 
hand 


(82) Enya = Volm — Vite 
say. On the other hand, it is easy to see that 


If now we put 


(8) and (9) evidently imply 
D; = V; = Voo<™ = 


where g;‘” is the jth elementary symmetric function of the quan- 


tities x°", (i=0, - - - , m—1). Then by (5) and (7) we have the 
formula 
(10) = (= — 


Returning now to (4), we see that 
= TT {fi + 


and therefore, making use of (10), we have finally the general- 
ized identity 


(i=0,---,m—1), 
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i(é). 


o{f™ 


(11) ( 


For 7 =0, this reduces to the identity (2). 


3. Necessity of the Criterion for Solvability. Returning to the 
congruence (1), we see, by repeated applications of (2), that, 
for m>k=deg P, 


k) 


Yall) =o (mod P). 


In the case m=k, it is evident that 


v(t) =P" (mod P); 


and clearly the congruence 


A (mod P) 


is solvable only for A =0, in which case it has p** solutions. We 
may therefore assume without any loss in generality that m<k. 

Making use of (11), we note that if (1) is solvable, then also 
the congruences 


Fa 
(12) — u;=A, j= 0,---,m-—1), 
i i 


of") 


are solvable. As for (12), we shall need the following theorem. 
If we put 


4, 
then the congruence 


—u=A (mod P) 


is solvable if and only if the product AP’ is congruent (mod P) toa 
polynomial of degree <k—1.* 
Applying this criterion to (12) we have the following theorem. 


THEOREM 1. If (1) is solvable (m<k), then necessarily each 
product 


* Duke Mathematical Journal, vol. 1 (1935), p. 166. 
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of") 
(13) ( ) AP’, =0,---,m-—1), 


is congruent (mod P) to a polynomial of degree <k—1. 


4. Sufficiency of the Criterion. Making a slight change in nota- 
tion, we rewrite (11) in the form 


(14) = — a;f(t), a;= ( ) ; 


we put 
(15) fi = (j = 0,---,m-— 1). 
Let us now consider the system of congruences in fo, « - - , fm—1: 
(16) =; (mod P). 


i=0 


The determinant B=|8;;| is easily calculated. By (14) we have 
the recursion 


m 
— = | Bj, = 1. 


Then 


B” =| =| (4,7 =0,---,m-—1), 


m m 
| |+ | | + a$;,n-1, 1 | 
1 m— 1 


| , 1 | 
Fy, 

* * Am—1 


(FiF2-- 


? 


and therefore 
: 
(17) B= ’ (c in GF(p")). 


0102°** Gm-1 FiF2 


= 
= 
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In particular B40 (mod P), and hence for fixed u; the system 
(16) has a unique set of solutions #;. It remains to show that, 
if the u; satisfy the congruence (12), then 


(18) fin , (¢=0,---,m—2). 


The proof follows closely the proof of (17). For simplicity we 
take only the case i =0 of (18); the general case may be treated 
in exactly the same way, but the notation is somewhat more 
complicated. We have 


Bto =| 4; A 
(19) 9 | i il i 1| (j = 0,---,m—1). 
Bi = | Bjo 
Then 
Bite = 87, 1 | 
m m | 
=/A a ~ a ;B j2, a j,m~1, 1 
2 m—1 | 
La 
= a jo, a jj, ;B j2, | 
Ln 
= — An 1 


Comparing with (17) and (19), we see at once that 4; =f)". We 
have therefore proved the following result. 


THEOREM 2. The criterion (13) furnishes a set of sufficient con- 
ditions that the congruence (1) be solvable. To carry out the solution 
we first solve the congruences (12) for u;, substitute these values in 
the right member of (16), and solve for t;=t". 


Note finally that the general solution of (1) is t+G, where t 
is a particular solution and G is any polynomial of degree <m; 
this follows from 


+ G) = + = 
5. Another Criterion. It is easily seen that the congruence 


(mod P) 


— 
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is solvable only when 
(20) A+AM+--- 4A", 0, 


In any event the sum (20) is congruent to a quantity in GF(p"); 
we denote this residue by p(A). Thus the congruence (12) is 
solvable only when 


{m—1)\ p” 


m—1 


Consider now the sum 


m—1 

j=0 
m—1 

i=0 
k-1 

j=0 

k-1 

i=0 


where for brevity we put 
A 


= M = 


Clearly (21) implies p,,(1/) =0. Conversely, it is easy to see that 
if pm =0, then (21) must hold. Indeed, since the sum (22) is con- 
gruent to a polynomial of degree m—1 (which is surely <), 
its vanishing entails the vanishing of its m coefficients; these are 
precisely the left members of (21). Thus we have the following 
theorem. 


THEOREM 3. The congruence 


= MF, (mod P) 


= 
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is solvable if and only if 


pm(M) = >> = 0. 
i=0 L; 


It will be remarked that in general the first criterion is more 
useful. However, the latter is of some interest in itself. Further- 
more, it suggests possible criteria for more general classes of con- 
gruences; I hope to develop the matter elsewhere. 


DvuKeE UNIVERSITY 


A DIFFERENTIAL EQUATION FOR APPELL 
POLYNOMIALS{ 


BY I. M. SHEFFER 


By a set of polynomials {P,(x)}, (n=0, 1, 2, - -- ), we shall 
mean an infinite sequence in which P,(x) is of degree exactlyf n. 
Corresponding to a given set {P,} there are infinitely many 
sequences of polynomials {L,(x)} (with L,(x) of degree not ex- 
ceeding ) and sequences of numbers {),} such that { P,} satis- 
fies the linear differential equation (usually of infinite order) 
with parameter :§ 


(1) L[y(x)] = La(x)y™(x) = dry(x), 


which for \=X, gives P,(x). In fact, suppose {P,} is given. 
Let {X,} be any sequence of numbers subject only to the condi- 
tion that X, is not identically zero in nm. Then a unique sequence 
{L,(x)} exists such that L[P,(x)]=\,P,(x), (n=0, 1,---), 
where not all the L,’s are identically zero, and where no L,(x) 
is of degree exceeding n. The polynomial L,(x) is readily ob- 
tained by recurrence from Lo, - - - , Ln-1. If we write 


(2) = Ino + + + Inne”, 


+ Presented to the Society, April 25, 1935. 

t For many purposes it suffices to have P,(x) of degree not exceeding n. 
Here, however, it is convenient to use the stricter condition. 

§ See Sheffer, American Journal of Mathematics, vol. 53 (1931), pp. 29-30, 
for a relation suggestive of (1). 
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where /,,, and other coefficients may be zero, we shall have 
(3) An = loo + + n(n 1)le2 + + 


If we start with equation (1), that is, with the coefficients 
{L,(x)}, and ask for what values of \ there exists a polynomial 
solution, we find as a necessary condition that \ have one of the 
values (3). Conversely, if \ has any value i, of (3), then there 
will correspond a polynomial solution P,,(x) which is unique to 
within a multiplicative constant, provided \..+\,, (mn). 
Moreover, P,,(x) will be of degree exactly , so that {P,} isa 
set. In the case in which \,, =A, for two (or more) distinct values 
of m and n, there may or may not be, for this common charac- 
teristic number, as many linearly independent polynomial solu- 
tions as the order of the characteristic number. 

DEFINITION. Let us suppose that equation (1) is given. If for 
each i, given by (3) there is a polynomial solution P,(x) of 
degree (exactly) m, we shall say that equation (1) is non-singular. 


COROLLARY. Equation (1) is non-singular tf \m¥#Xn, (mM¥n). 


If equation (1) is non-singular, but \,,=A, for some m#n, 
then if P,, and P,, are corresponding solutions of degree m and n, 
the polynomial aP,,+5P, is also a solution, for every choice of 
a and b. In this case the equation does not define a uniquet set 
of polynomials { P,,}. 

Because of the property of equation (1) that amy set { P(x) } 
is a set of solutions, (for a suitably chosen sequence {L.(x) } ); 
it seems appropriate to refer to (1) as a universal type equation 
for sets of polynomials. This comprehensive character of a type 
equation suggests that it may prove useful in the study of prop- 
erties of sets of polynomials. We may propose the following pro- 
gram. Given a set of polynomials, or a class of sets; to charac- 
terize the type equations that they satisfy, and, in particular, to 
determine if (and when) they correspond to a type equation of 
finite order. 

The present note investigates these questions for the class of 
Appell sets of polynomials. 

The set {P,,} is an Appell set if we have the relations 


+ We ignore the obvious arbitrary multiplicative constant attached to each 
P,,(x). 


| 
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dP,,(x) 
(4) = P,_1(x), (n = 1,2,---). 
dx 


The following are well known characterizations of Appell sets. 
(i) A necessary and sufficient condition that {P,} be an Appell 
set is that a formal power series 


(5) A(t) ~ ant” 
0 

exist such that 

(6) >> P,(x)t*. 
0 


We shall refer to A(t) as the generating function (or power series) 
for { P,(x)}. 
(ii) Set { P,.} isan A ppell set if and only if there exists a sequence 
of numbers {a,} such thatt 


(n = 0,1,---). 


THEOREM 1. Let [P,3 be any Appell set, with generating func- 
tion A(t). Define the (formal) power series 


(8) Bit) ~ > bt” 
by 

A'(t) 
(9) B(t) = A 


Then {P,,} satisfies the equation of type (1): 


+ In both (6) and (7), in order that {P,} be a true set, that is, for P, to 
be of degree exactly n, it is necessary (and sufficient) that a9~0. There is no 
great difficulty, however, if aa= --- =a,=0, For then P,(x), as given 
by (6) and (7), is identically zero for OS" Sk, after which we get P,(x) of 
degree n —(k+1); and if we define B(t) by A =t**1B(t), then bo +0, and the 
Appell set for B(#) is the same set as for A (é) after the first (+1) P’s (which 
are identically zero) for A(t) have been omitted. 
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L[P,(x) | 
(10 


(bo + x)Px (x) + (x) + (x) +--- 
= nP,(x), (n = 0,1,---). 


_ Conversely, if the set {P,} satisfies the equation (10), where 
bo, bi, - - - are constants, then (to within multiplicative constants) 
{P,} is an Appell set, and its generating function A(t) is defined 
by (9). 


From (3) we see that 4,=m for equation (10). Therefore 
Am*~An, (mn), so that toeach i, corresponds an essentially 
unique polynomial, and it is of degree n. Suppose the set { P,} 
satisfies (10). Let Po{x) =ao, and choose the multiplicative con- 
stants so that in P,(x) the coefficient of x” is ao/n!. It is 
verified that P/ (x) = Po(x). Assuming that P/ (x) =P for 
k=1, 2,---,mn-—1, we shall prove it for k=n. Differentiating 
(10), we find that L[P,! (x)]=(n—1)P,/ (x). Hence, since L is 
non-singular, we have P,! (x)=c,P,:(x). Comparing coeffi- 
cients of x"~', we find that c,=1. Hence {P,} is an Appell set. 
Let its generating function be A()~ vant”, where Py as 
already stated. Then P,,(x) has the representation (7). If we sub- 
stitute this in (10) and equate coefficients of lke powers of x, 
we get the equations 


+ + + Oo = (8 — 
VG =0,1,---,m). 


But these relations are completely equivalent to relation (9), 
as is easily verified. The converse part of Theorem 1 is thus es- 
tablished. 

Now let { P,,} be any Appell set, and let A (¢) be its generating 
function. Let us define B(t) by (9), so that the coefficients of 
(10) are defined. Let {Q,} be the Appell set defined by (10), 
of whose existence (and uniqueness if we choose Qo(x) =a) we 
have just given proof. Then the generating function A*(¢) for 
{Q,} is defined, by the part already proved, by the equa- 
tion B(t) =A*’(t)/A*(t). Since A(0)=A*(0)=ao, therefore 
A*(t)=A(t). Hence Q,(x)=P,(x), (n=0,1,---), so that 
{P,,} satisfies (10). This completes the proof. 


COROLLARY. A necessary and sufficient condition that a set 
ipa be an Appell set is that it satisfy the recurrence relations 
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(11) nP,,(x) (bo + Pp-1(x) + bP ,-2(x) b,-1Po(x), 
(s = 0,1,---). 


For if (P.} is an Appell set, it satisfies (10) and therefore (11). 
Now suppose {P,} is a set satisfying (11). We verify at once 
that =P». Suppose P? = Pi, (k=1, 2, - - - ,m—1); to prove 
that P,’ =P,_,. If we differentiate (11) and use our induction 
assumption, we get wP,{ =nP,_1, so that P,’ = P,_1, and the set 
{P,} is an Appell set. ; 

Relations (10) and (11) give us two further characterizations 
of Appell sets. But they are not the only ones possible. For ex- 
ample, if we want an equation (for Appell sets) of type (1) with 
\, =u(n—1), then as was proved in Theorem 1, we obtain the 
equation 


(12) = n(n — 1)P,(x), 
where the operator Ly) is defined by the equation 


Lely] = Do (x) 
with 


(13) 2) ~ DS = + 2xA’ + x?A], 
n=0 


where A is the generating function for 
More generally, if { P,} is an Appell set defined by A(t), then 


(14) Luy[P.(«)] = n(n —1)--- (n — k +1)P,(x), 


where 
k 


t k 
1) ~ >. = + 2A 
(15) n=0 


Now suppose that we wish X, to be a polynomial in . We can 
suppose thatf A» =0. Then we have 


An = + n(n — + 


+ For this merely serves to alter the (constant) value of Lo(x), and does not 
affect the polynomials L,,(x) of index n>0. 


(16) 


+n(n—1)---(#—k+ 
| 
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where J, - - - ,/,, are constants (that is, independent of m). Ex- 
actly as in Theorem 1, we obtain the following result: 


THEOREM 2. Let {P,} be any Appell set and A(t) its generating 
function. Let {d,} be given by (16) for any k=1. Then { P,,} satis- 
fies the equation 


(17) L[P.(2)] = = rnPa(2), 
r=0 

where L is defined by 

(18) L(x; t)~ L,(x)t* 


= (x; t) + t) + - + 2), 


the functions Li(x; t) being given by (15). Conversely, given (17) 
with L defined by (18) for any function A(t), then there exists an 
Appell set of polynomialst {P,} satisfying (17), and its generat- 
ing function is precisely A(t). 


We proceed to characterize the equations of type (1) that de- 
fine one and the same Appell set. 


TueoreM 3. Let {P,,} be an Appell set, and A(t) its generating 
function. A necessary and sufficient condition that the set {P,} be 
a solution of an equation of type (1) (with Xx =0, which is no re- 
striction) is that the operator L Ly] (x) be definedt 
by 


(19) 1) ~ ~ (x; 2), 
n=0 


where Lyi(x; t) is given by (15), and the l,,’s are constants. The 
characteristic numbers {d,} are given by (3). 


If Xm An, (mn), then this Appell set gives al] the polynomial solutions 
of (17) (to within multiplicative constants). But if mand m (distinct) exist such 
that A» =An, this Appell set does not include all the polynomial solutions of 
(17). 

t Although the last series in (19) may not converge, yet it serves to deter- 
mine a unique sequence of polynomials L;(x), L2(x), - - - which are the coeffi- 
cients of equation (1). For L,x)(x; t), expressed as a power series in ¢, begins with 
a term in ¢*. Hence L,(x), which is the coefficient of ¢” in L(x; t), is given by a 
finite sum, namely, the coefficient of in the sum (x; 2). 


| 
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First we shall show that for any sequence of numbers {ix} 
the Appell set { P,,} satisfies equation (1) with L defined by (19). 
Let 


[n] = te), 
(%3t) = DL, (x), 
r=0 


k=0 


where 1 is any positive integer. The remark of the last footnote 


shows that =L,(x), (r=0, 1, - - - , 2). Now, by Theorem 
2, we have 
[P,(x)] = {rl + — + + Pr(x) 
= ),P,(x), (ry =0,1,---). 
Hence 


L[P,(x)] [P,(x) ] \-P,(x), (r 0, 1, n); 


and inasmuch as m can be chosen arbitrarily, it follows that 
L[P.(x)]=\,P,(x) for all r. That is, { P,,} satisfies equation (1). 

Conversely, suppose the Appell set {P,} satisfies equation 
(1) (with \» =0) ; we shall show that ZL is defined by (19) for suit- 
ably chosen numbers {/..}. Define {lx} by (3), and with this 
choice of {J:.}, define an operator L* by (19): 


L*(x;t)~ 2), 
k=1 
where the function A(t) of (15) is the generating function for 
{P,,}. From the half of the theorem already established, { P,} 
satisfies the equation L*[y]=)*y, the set {d,*} being given by 
(3). But {X,} also is given by (3), so that \,* =,. Subtracting 
the two equations, we find 


(L — L*)[P,.] = 0, 0,1,---). 
If we set n=0, 1, - - - successively, we find that 
Lo(x)—Le'(x) =Li(x)-Li(x)= --- =L,(x)—L,*(x)= --- =0. 


That is, L*=L. This proves the converse. 

We next determine the condition that an Appell set satisfy 
a finite order equation of type (1). Two well known examples are 
the polynomial sets { (x —a) and { H,(x) }, the latter being 


| 

a 
i 
} 
| 


1935-] APPELL POLYNOMIALS 921 


the Hermite polynomials. They satisfy the following equations, 
respectively: 


L[y(x)] = — a)y'(x) = dy(z), = 1); 
L[y(x)] = 2ay'(x) — = Ay(z), On = 222). 


The corresponding generating functions are respectively 
A(t) and A(t) =e-*!/4, 


THEOREM 4. A necessary and sufficient condition that an A ppell 
set {P,, } , with generating function A(t), satisfy a finite order equa- 
tion of type (1) is that A(t) have the form 


(20) A(t) = 
where Q(t) is a polynomial. 


Suppose that A(#) has the form (20), where Q is of degree k. 
Then from (9), B(t)=Q’(t), a polynomial of degree k—1, so 
that B(t)=bo+bit+ - - - The universal equation (10) 
for Appell polynomial sets therefore reduces to 


L[Px(x)] = (bo (x) + +--+ + (2) 
= nP,(x), 


and this is an equation of finite order (namely, order k). 

Conversely, suppose the Appell set {P..} satisfies a finite 
order equation of type (1); we shall show that its generating 
function A(t) has the form (20). Let the equation be 


L[Pa(x)] = Lo(x)Pa(x) + (x) = 
= 0,1,---). 


As already pointed out, we may suppose that \»=0 (so that 
Lo(x) =0). If we write L,(x) in the form (2), then X, is given by 
(3), but since L,(x)=0, (n>k), \, reduces to the form (16). 
Let Lyy(x; t), (k=1, 2,--- ), be the functions of (15), and de- 
fine the operator L* by L*(x; #) a Li (x; t), where A(t) 
is the generating function for {P,}. Then {P,} satisfies (by 
Theorem 2) the equation L*[P,]=,*P,, where \,* is deter- 
mined by (16). Since (16) was used to determine X,, it follows 


| 

| 

Ca... ... 
i 

| 

i 

| 

| 
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that A,*=X, for all m. Hence (as was argued before) we have 
L*=L. That is, the given equation for {P,} is defined by 


k 
(a) L(x; #) = 


But to begin with we had 


k 
(b) L(x; t) = L,(x)e. 
r=1 

Hence 

k k 
(c) t) = L,(x)t". 

r=1 r=1 
Let 

A(t) 

(d) (r= 0,1,---,k). 


A(t) 


Expanding the left member of (c) in a power series in ¢, and 
observing (from the right member) that powers of ¢ higher than 
the kth degree cannot occur, we obtain the equations 


+ + 2%01,n-1) °° 


k 
+ lex (cries + + 


(e) 
| 
(k — 1)! 


= 0, (n = k). 


These relations are to hold identically in x. Of the numbers 
1,---, &, let s be the largest for which /,,~0. It then follows 
from (e) that 


(f) @1,n-s+1 = 0, (n = k). 
That is, 


— = Bit) = t 
A) ( ) + + + 41% 


__ 
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so that (20) holds, where Q(¢) is a polynomial of degree k—s+1. 
This completes the proof. 

If the Appell set {P,.} satisfies an equation of type (1) of 
finite order, then there are infinitely many other equations (of 
type (1)) of finite order for which the set is a solution. For ex- 
ample, from follows L?[P,,]=,2 P», Pa, 

- + +, This raises the question : What is the minimum order of all 
equations satisfied by { P,}? The answer is given by the follow- 
ing statement. 


THEOREM 5. Let the A ppell set {P.} satisfy a finite order equa- 
tion of type (1), so that its generating function A(t) 1s given by (20). 
Then the minimum order of all equations (of type (1)) satisfied by 
{P,} is precisely the degree} of the polynomial Q(t). 


Let Q(t) be of degree m. That the minimum order k cannot ex- 
ceed m is shown by the proof of the first half of Theorem 4, so 
that k<m. Lett L[P,]=\,P, be the equation of minimum 
order k. From the proof of the second part of Theorem 4, we 
see that there is an integer s for which 1<s<k, and such that 
Q(t) is of degree k—s+1. That is, m=k—s+1. But km, so 
that k<&k—s+41. This requires that s=1, from which follows 
the desired conclusion that k=m. 


PENNSYLVANIA STATE COLLEGE 


+ The single exception is when Q(f) is of degree zero, in which case A(t) isa 
constant and P,,(x) =cx"/n!. This set satisfies a first order equation. 
t As usual it is no restriction to suppose, as we do, that A»=0. 
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ON SOME EXTREMAL PROPERTIES OF 
TRIGONOMETRIC POLYNOMIALS 
WITH REAL ROOTS 


BY J. GERONIMUS 


1. Introduction. L. Fejér [1],f O. Sz4sz [2], [3], [4] and E. 
v. Egervary [5] have found many interesting extremal prop- 
erties of non-negative trigonometric polynomials. In particular, 
Szasz [3] has found that for every non-negative trigonometric 
polynomial of order <n with real coefficients, 


(1)G,@)=14+RQ yee, (ye = oe + = 1,2,--- , m), 
k=1 
the inequality 
us 
(2) | < 2 cos ———_-, (k = 1,2,---,#), 


k 
is valid.f 


The object of this note is to find the minimum of the modulus 
of the first coefficient y,, supposing that all roots of G,(@) are 
real. The first problem of this kind has been considered by 
Blumenthal [6]; we shall return in §4 to his problem and its 
generalization. 


2. Equality of Roots of G,*(@) for Problem 1. Consider the fol- 
lowing problem. 


PROBLEM 1. Find the minimum of the modulus of the first coeffi- 
cient Y, of a non-negative trigonometric polynomial 


G,(0) =1+R> 


k=1 


of order n with real roots. 


{7 Numbers in brackets refer to the Bibliography at the end. 
t Rz means real part of z; [a] means the greatest integer <a. 
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In order to solve this problem we shal! prove the following 
simple lemma. 


Lemna 1. All roots of the polynomial G,*(0) for which the mini- 
mum in Problem 1 1s attained must be equal. 


Consider a non-negative trigonometric polynomial t 


A 


0 — = 
(3) G,(0) = 4 sint —— sin? — F,-2(6) = Frei, 
k=0 


where F,_2(@) is a non-negative trigonometric polynomial of 
order 2 —2 with real roots, 


n—2 n—2 

(4) Fy-2(0) = RY Vee = | cos (k0 — ax), 
k=0 k=0 

where a; =arg y*, 


On putting a= (0:+62)/2, 5=(6:—62)/2, we see easily that 


1 
Yo = + cos 28) — | y#*| cos (a1 — a) cos 
(5) : 
+ | cos (az — 2a); |v] = 
We see that |7,| does not depend on a, nor on 6; on the other 
hand yo is maximal for 6=0 if cos (a,—a) $0, or for 6=7 if 
cos (a:—a) 20. In both cases the minimal value of |-y,,| under 


condition yo =1 corresponds to 6=0 or 6=7; therefore and 6, 
coincide. f 


3. Polynomials for which y, has Extremal Values. It follows 
from this lemma that G,*(6) is 


(6) G*(0) = C[1 + cos (9 + a)]", 


a being an arbitrary real argument; it may be written thus:§ 


(7) G.*(0) = 


1 n 
Conn + Conn—z cos + : 


k=1 


Qn-l 


{ It isclear that all real roots of a non-negative trigonometric polynomial 
are of even multiplicity. 

t 6, and 6,+27 are not considered as different. 

§ See [6], p. 392. 


= 
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For this polynomial we have 
Qn Qn-1 
whence we find the ratio 


2 
Yo Cao 


(8) 


We have proved the following theorem. 


THEOREM 1. If G,,(8) is a non-negative trigonometric polynomial, 
G0) =1+ RY 
k=1 


of order n with real roots, then 


2 


©) Canin 


< |y,| <1; 


the maximum is attained for the polynomial} 
(10) Ginax(@) = 1 > cos a), 
and the minimum for the polynomial 


n 


(11) Gmin(0) = }1 + cos (0+ a) " 


2n,n 
a being an arbitrary real argument. 


4. The Generalized Extremal Problem. Consider now the fol- 
lowing extremal problem. 


PROBLEM 2. Find the minimum of the ratio 
+ 
(12) 
A? + (AP + Be) 


k=1 


t See [1], [2]. 
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where 

gm(9) = Ao + Ai + B, + --- + A, cos m0 + B, sin 
is a trigonometric polynomial of order m with real roots, and d is an 
arbitrary non-negative number. 


The above mentioned problem of Blumenthal [6] corresponds 
to \=1. It is easy to see that for \=2 Problem 2 is a particular 
case of the Problem 1. Indeed we see that 


(13) ge (0) = = RYO 


k=0 


is a non-negative trigonometric polynomial of order »=2m, 
while 


i> A,2 + B,2 
(14) =A? +— (42 + B2), | = ————; 

therefore we have for \==2 

+ B,2 2 2 
(15) 

(42+ Be) 
k=1 


To solve our problem for all X20 we,shall put it in the following 
form. 


PROBLEM 2’. Find the maximum of the expression 


1 2r 1 2 
(16) LG) = + f IG,,(0) | ah, 


(e 1), 
where G,(0) 1s a non-negative trigonometric polynomial 
G,(0) = RY = 1), 
k=0 
of order n=2m with real roots. 


5. Equality of Roots of G,*(0) for Problem 2’. We shall prove the 
following lemma. 


= 
= 
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LEMMA 2. All roots of the polynomial G,*(0) for which the maxi- 
mum in Problem 2’ is attained must be equal. 
Put 
A, 0 — 
G,,(0) = 2 sin sin Fn—1(0 
(17) IG,.(6) ] 1(8) 
[cos 6 — cos (@ — a) |Fn_1(8), 


where a= (0,+62)/2, 6=(0:—62)/2, and Fn—i(@) is a non-nega- 
tive trigonometric polynomial of order m—1, 


(18) = RY 


with real roots. Thus we get 
1 ae 1 

(19) —{ IG,,(0) = co cos 6 — —| cos (8; — a), 
0 2 


where =arg cz, (k=0, 1, --- , m—1), and Further let 


2m— 


(20) Fn-1(0) =R 


then we obtain 


1 al 1 
G,,(0)d@ = +— cos 28) 
0 2 


(21) 
1 
| cos — a) cos6 + | cos (6s — 2a), 


2 
=0 


where 6*=arg ci*, (k=0, 1,---, 2m—2), and B*=0. Using 
(19) and (21) we have 


(22) LG) = A cos 26+ Bcos6+C, 
where 
A= (co* + 
B = —|c#| cos (8* — a) — c1| cos — 


(23) 1 
C= = c#| cos (B# — 2a) 


1 1 
+ Jel (B1 — a) + + 
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It is important to point out that |y,| =(1/4)| c2m-2| does not de- 
pend on a, nor on 6. Since we have 


1 m—1 
(24) cr = +—Dla|?, 
2 
it is clear that for «= —1 we have 
1 
(25) 
2 4 
Therefore L(G,) is maximal for 6=0 if B20, and for 6=7 if 


B<0; in both cases 6; and @, coincide, which proves our lemma. 


6. Polynomials having the Extremal Property. We see that the 
polynomial G,*(@) is 


G.*(0) 


2™-1[1 + cos (@ + a) |” 
26 1 ; 
( ) = Cann + cos + a), 


k=1 


and we have for it 


(27) LG,*) 


1 


Thus we have proved the following theorem. 


THEOREM 2. If G,,(0) is a non-negative trigonometric polynomial 
of order n=2m, 


G,,(0) R>, | | = 
k=0 


with real roots, then 


1 Qr 1 Qr 2 
Gu(0)40-+ — G.(0) ] nao 


(28) ‘ 
Ss 2 + €(Cn,nj2)”) (e 1); 


the maximum is attained for the polynomial 
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(29) = 2"-1{1 + cos (6 + «)}”, 
a being an arbitrary real argument. 
This result may also be stated as the following theorem. 
THEOREM 2’. If gm(@) is a trigonometric polynomial of order m, 
gm(0) = Ag+ Aicos@+ --- +A, cosmé+ B,, sin 


with real roots, then 


A,2 + B,2 2 
(30) = » (A2 0); 
Me + (AZ + BE) Cam,2m + Cem.n)? 
k=1 


this minimum 1s attained for the polynomial 
(31) = C{1+ cos (0+ a)}™. 
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LINEAR CONNECTIONS OF NORMAL SPACE TO A 
VARIETY IN EUCLIDEAN SPACE* 


BY C. B. TOMPKINS 


1. Introduction. This paper deals with an extension of the 
Gauss formulas of a surface imbedded in ordinary space to ap- 
ply to an m-dimensional variety imbedded in an n-dimensional 
euclidean space. The Gauss and Codazzi relations are extended 
to give integrability conditions in terms of the Christoffel sym- 
bols of the second kind of the variety, a set of tensors corre- 
sponding to the coefficients of the second fundamental form of a 
surface in three space and a set of non-covariant quantities: the 
connection coefficients of the normal space. 

Several authorst have considered generalizations of the 
Frenet or the Gauss formulas to apply to varieties lying both 
in euclidean and in more general spaces, but the nature of some 
of the coefficients which may appear in these formulas seems to 
have escaped serious study. 

In the considerations here there will be encountered an im- 
bedding euclidean space of m dimensions, a variety of m dimen- 
sions lying in it and at each point of this variety an m-dimen- 
sional tangent space and an (m—m)-dimensional normal space. 
A vector in the imbedding space will be denoted simply by a 
letter, and all indices running from 1 to m will be suppressed. 
Latin indices lying between a and k will have the range from 1 
to m and Latin letters from p through z when used as indices 
will have the range from 1 to (n—m). Corresponding Greek let- 
ters will be used as summation indices. : 

The work in this paper differs from that in most of the previ- 


* Presented to the Society, April 19, 1935. 

T See, for instance, Voss, Mathematische Annalen, vol. 16 (1880), p. 129; 
H. Weyl, Mathematische Zeitschrift, vol. 12 (1922), p. 162; Schouten and van 
Kampen, Mathematische Annalen, vol. 105 (1931), p. 144; E. Bortolotti, 
Rendiconti del R. Istituto Lombardo di Scienze e Lettere, (2), vol. 64 (1931), 
p. 441; E. H. Cutler, Transactions of this Society, vol. 33 (1931), p. 832; C. E. 
Weatherburn, Reports of the Australian and New Zealand Association for the 
Advancement of Science, vol. 21 (1933), p. 12; Duschek and Mayer, Lehrbuch 
der Differentialgeometrie, vol. 2. 
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ously published papers in the method of selecting coordinates 
in the normal space; in this paper a set of perfectly arbitrary 
mutually orthogonal unitary vectors is set up as a coordinate 
system in the normal space at every point of the variety, while 
in most of the other papers the normal space is divided into 
several subspaces (the first containing all components of the 
first derivatives of the tangent vectors lying in the normal space, 
the second determined by additional components of the second 
derivatives, and so on) and in each of these subspaces an arbi- 
trary coordinate system is introduced. Weyl derived the same 
Gauss formulas and integrability conditions we shall use, but he 
seems not to have noticed the properties of the coefficients which 
we shall call N%&. 


2. The Gauss Formulas and Integrability Conditions. A variety 
of m dimensions may be represented by a set of equations de- 
pending on m parameters, x=x(u*). In general there is a 
uniquely determined flat m-space tangent to this variety at 
every point; this is the linear space depending on the vectors 
p:=0x/du'; it is m-dimensional and uniquely determined if the 
vectors p; are linearly independent. Also at every point there is 
a flat space which is absolutely perpendicular to the tangent 
space. If the tangent space is actually m-dimensional, the nor- 
mal space is (n—m)-dimensional and contains sets of (n—m) 
unit mutually orthogonal vectors. We may choose any such set 
and name it ¢,. Now, any vector connected with a point of the 
variety may be written as a linear combination of the vectors 
p; and ¢t, associated with that point. In particular we can write 


Op: 


= + Lid, 


(1) 


where the coefficients Ti; are the Christoffel symbols of the sec- 
ond kind connected with the variety and the choice of parame- 
ters. 

We write also the equations for the partial derivatives of the 
vectors f,, 


aly 
(2) Li” pa + N pity, 
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where it is easily shown that L, =g,,.Li”. Because the vectors 
t, are mutually orthogonal unitary vectors, the coefficients N{% 
must be anti-symmetric in the indices p and gq. 

The simplest method of arriving at integrability conditions 
seems to be that of variation of product integrals as developed 
by Rainich and Vaughan.* It is well known that we can take 
any set of equations of the form of (1) and (2) and by applying 
product integrals over some path arrive at a set of values for 
the vectors p; and ¢, along the path. However, for the variety 
to be a true variety it is necessary and sufficient that the vectors 
obtained for any point be the same no matter what path is 
used in deriving them; that is, the variation of the product in- 
tegral used in their derivation must be zero. If we denote the 
matrices of the coefficients of the ; and ¢, in the expression for 
the partial derivatives with respect to u* in (1) and (2) by Ai, a 
necessary and sufficient condition for the variation to vanish is 

0A; OA, 


A;Ay — t+ = 0 
Out 


If we denote the matrices which make up A, in the obvious 
manner by T;, Li, Li’, and N;, three relations analogous to the 
Gauss and Codazzi relations follow, conditions necessary and 
sufficient for the set of equations (1) and (2) to represent a true 
variety: 


or; or; 
Ry; = — = = 
ou* out 
} ON; ON: 
(3) = — —-+ — NiN; = — Lili, 
OL; OL; 


+ + LiN: LN; 0. 


ou* oui 


Having given the quantities [',, Li, and N; satisfying these con- 
ditions and having given also a set of values for the components 
of the fundamental metric tensor at some point on the variety, 
the variety may be completely reconstructed except for position 
by means of a product integration followed by a Riemann inte- 
gration. 


* Abstract, this Bulletin, vol. 40 (1934), p. 233. 
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3. Laws of Transformation. All the coefficients of the equa- 
tions (1) and (2) are invariant under transformations of the co- 
ordinate system of the imbedding euclidean space. Under trans- 
formations of the system of parameters u*, the coefficients Tf, 
transform under the well known law for the transformation of 
connection coefficients; under these transformations, the coeffi- 
cients Li. transform as second rank doubly covariant tensors 
with indices 7 and k, while the similar quantities with one index 
raised must obviously transform as second rank mixed tensors. 
These transformations are all obvious from the form of equation 
(1). From equation (2) it is also evident that under arbitrary 
changes in parameterization of the variety the coefficients Ni, 
transform as first rank covariant tensors with index k. 

In the normal space we allow rotations; more general trans- 
formations could be introduced with slight changes in the for- 
mulas already developed, but the introduction of such additional 
generality would simply necessitate the distinction between co- 
variance and contravariance without rendering any useful serv- 
ice. We consider transformations of the type /=At, where A is 
an orthogonal matrix. Applying this to the quantities of equation 
(1), we get the law of transformation for the coefficients Li,, 

Liu = 
By continuing to equations (2), we find the law of transforma- 
tion of the other coefficients, 


oA? 


ou* 


= AN,AT —A 


This law of transformation is formally exactly the same as the 
law which must be obeyed by coefficients of linear connection. 
This permits the statement of the following theorem concerning 
differentiation of complete tensors.* 


THEOREM. The operation of covariant differentiation on a com- 
plete tensor,in which operation the differentiation is with respect to 
one of the parameters u*, and in which the quantities Nj, are used 


* That is, a set of quantities which are components of tensors in all the 
coordinate systems we are considering; the process we describe has been called 
complete differentiation (see E. H. Cutler, loc. cit.). 
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as connection coefficients for indices applying to the normal space, 
and in which the ordinary Christoffel symbols, rs, are used as 
connection coefficients for indices appertaining to the tangent 
space, yields a complete tensor. 


The expression S}; of the second of the integrability condi- 
tions (3) is simply the Riemann tensor of the connection coeffi- 
cients of the normal space. The analogy between the first and 
the second of these conditions is striking. The third of these in- 
tegrability conditions is rewritten under the definition of covari- 
ant differentiation just stated as 


Ads; = 0. 


The anti-symmetry of the connection coefficients of the normal 
space reiterates the equivalence of covariance and contravari- 
ance. 


4. Vanishing of the Riemann Tensors. It may be that there 
exists a choice of the vectors ¢, such that the coefficients N%, 
vanish simultaneously; from the law of transformation of these 
coefficients, it is evident that a necessary and sufficient condition 
for the existence of such coordinates is that the coefficients have 
the form N; = —A(0A7/du*); or, since the matrix A is orthogo- 
nal, 0A /du* = N;A. This equation is of the type whose solution 
for the matrix A may be expressed by product integrals; we can 
find a solution over any particular path, but the solution must 
not depend on the path; that is, the variation of the product in- 
tegral must vanish. A necessary and sufficient condition for this 
is the vanishing of the Riemann tensor of the normal connection 
coefficients, Si.;=0. The solution we obtain by means of such 
an integration for A is an orthogonal matrix, for it can be shown 
that the product integral of any anti-symmetric matrix is or- 
thogonal if the matrix used as an arbitrary constant of integra- 
tion is orthogonal. Hence, the vanishing of the Riemann tensor 
Shi 1S a necessary and sufficient condition for the existence of co- 
ordinate systems for which all the N}, vanish. 

The implications of the vanishing of this tensor are readily 
seen from the second of the integrability conditions (3); if we 
choose the parameters u‘ in such a way that at a point in which 
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we are interested the p; are unit mutually orthogonal vectors, 
then L;’ = L{;, and the integrability condition becomes 


Lig 
Lis Lie 


Now rotate the vectors so that one of the matrices, say L’., be- 
comes a diagonal matrix, and the equation becomes 


— Lin) = 0. 


This means that the vanishing of the tensor Sh; implies that at 
every point on the variety there exists a set of mutually orthogonal 
vectors in the tangent space of such a nature that all the principal 
directions of all the Dupin indicating quadrics lie along these vec- 
tors. This is evident when we notice that the equation implies 
either that the other matrices L*‘, are all diagonal when L’. is 
or that Li;=Ly. If this last condition is true, the matrix L’, 
has no principal directions with respect to the indices 7 and k, 
and we may rotate the vectors p; and #; in their plane until the 
matrix L*, becomes a diagonal matrix without destroying the 
diagonal property of L”.. 

In the case of a surface in a euclidean four-space we can get 
more definite results. The vanishing of the tensor Sj; implies 
that the indicatrix of Wilson and Moore* degenerates into a 
straight line. Indeed, if we choose the coordinate axis y along 
the vector ¢, and the axis z along é2, the indicating conic may be 
written parametrically 


y= Li: cos? + 2L sin cos + Liz sin? 6, 


cos? + sin 9 cos 6 + sin? @. 


Because of integrability conditions and the vanishing of the 


tensor S%,;, these reduce immediately to the equation of the line 
1 2 
— y—Lu -0 
1 2 
Li2 


UNIVERSITY OF MICHIGAN 


* Wilson and Moore, Proceedings of the American Academy of Arts and 
Sciences, vol. 52 (1916), p. 324. 
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